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INTRODUCTION 


As  previously  discussed  by  Kummler,  Fisher,  and  Boynton 
(1973) , hydvazine  and  substituted  amines  form  an  important 
class  of  fuels.  When  the  products  of  these  fuels  interact 
with  a low  pressure  environment  containing  oxygen  atoms,  mea- 
surable radiation  from  the  ultraviolet  through  the  infrared 
spectral  regions  occurs.  This  radiation  occurs  both  indirectly 
by  energy  transfer  to  infrared  active  products  of  the  substituted 
amine  combustion  (CC>2,  H20,  CO,  etc.)  and  directly  by  strongly 
exothermic  chemiluminescent  reactions  initiated  by  0 atom  attack. 
Most  of  the  energy  released  by  such  reactions  is  expected  to 
form  vibrationally  excited  reaction  products  and  is  therefore 
observable  as  infrared  radiation.  Minor  radiation  in  terms  of 
available  energy,  but  not  in  terms  of  observable  signals,  is 
expected  to  be  observable  in  the  visible  and  vacuum  ultraviolet 
spectral  regions.  Either  wavelength  region  provides  qualita- 
tive diagnostic  information  regarding  the  mechanism  of  reaction, 
but  the  value  of  the  diagnostic  depends  upon  our  quantitative 
knowledge  and  understanding  of  the  radiation  mechanism. 

Predictions  of  the  radiation  associated  with  chemilumines- 
cent processes  requires  knowledge  of  the  concentration  of  reac- 
tive species,  rate  constants  for  the  important  reactions,  and 
the  fraction  of  the  energy  liberated  by  the  reaction  which  ends 
up  in  a given  internal  mode.  Although  little  is  known  directly 
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about  unreacted  fuel  concentrations  in  ei  vironments  of  interest, 
there  is  some  evidence  that  levels  of  a few  percent  may  exist 
in  an  excess  of  atomic  oxygen. 

This  report  summarizes  the  experiments  conducted  in  the 
Environmental  Kinetics  Laboratory  of  Wayne  State  University 
which  demonstrate  the  chemiluminescent  processes  occurring  in 
these  systems.  The  systems  studied  include: 

0 + N~H . -*■  H o0  + HNNH 
2 4 2 


0 + NH. 


nh2  + OH 


0 + NR3  (where  R is  a hydrocarbon  radical) 
H + N02  -*•  NO  + OH  . 


products 


As  illustrated  in  Figure  1 (Kummler,  et  al,  1973)  all  of  these 
reactions  should  be  prominent  in  the  lean  fuel-oxygen  atom 
systems.  As  indicated  in  Table  I,  a preliminary  model  has  been 
formed  and  has  been  tested  in  simplified  form  (Kummler,  e t al , 
1973)  with  reasonable  assumptions  regarding  the  energy  distri- 
butions. In  this  work  we  attempt  to  support  those  assumptions 
by  experimental  work  on  individual  systems  in  the  context  of 
Figure  1 and  Table  I . 


HCH-NO 


FIGURK  1.  Simplified  O atom  oxidation  mechanism  for  MMH  (UOMH  and  n2H.  are 
similar)  illustrating  the  stable  species  which  can  form  in  this 
system  as  intermod iates . The  dotted  lines  represent  reactions 
which  ate  probaole,  but  which  have  not  yet  been  experimentally 
observed. 
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EXPERIMENTAL  APPARATUS  AND  PROCEDURES 

The  experimental  system,  illustrated  in  Figure  2,  has  been 
described  previously  in  related  publications  (Krieger,  Malki 
and  Kummler , 1972;  and  Kuirmler,  et:  al_,  1973)  . All  experiments 
were  performed  on  that  discharge  flow  tube  system.  Flow  rates 
up  to  80  m/sec  could  be  obtained.  Detection  systems  include  a 
quarter  meter  Jarrell  Ash  (M)  monochrometer  with  gratings  blazed 
at  3000  A,  6000  A,  2.1y,  and  5y.  An  RCA  C31025C  photomultiplier 
was  used  to  view  spectra  from  2000  A to  9000  A,  a Santa  Barbara 
ITO  Pbs  detector  from  1 - 3y,  a LN2  Au-Ge  detector  from  3-6y, 
and  a Santa  Barbara  Au-Ge  liquid  He  cooled  detector  was  used 

/ 

from  6 - 25y  (P) . A Raytheon  PGM  10  x 2 pulsed  diathermy  unit 
provided  a pulsed  discharge  at  30-10000  Hz  for  infrared  de- 
tection with  a PAR  FIR-8  phase  sensitive  lock  in  amplifier. 
Visible  detection  employed  a Keithley  602  electrometer  (E) . 
Honeywell  Electronic  194  dual  channel  recorders  (R)  were  used. 
Pressure  in  the  reactor  was  monitored  using  an  MKS  baratron  77H 
capacitance  manometer  with  a 10  torr  head.  C>2,  H2,  or  N2  were 
discharged  using  the  2450  MHz  diathermy  and  either  Evensen  or 
Broida  cavities.  An  Extranuclear  3/4"  quadruple  mass  spectro- 
meter was  available  to  the  system.  A Brad  Thompson  774  electron 
beam  welder  was  used  for  electron  excitation  studies.  Oxygen 
atoms  were  monitored  by  N02  titration.  Absolute  actinometry  was 
obtained  from  4000  A to  3y  using  the  NO  + 0 glow.  Gas  flows  were 


ig.2 -Schematic  of  Discharge  Flow  Tube  Reactor 

and  Detection  System 


monitored  by  rotometers  and  calibrated  capillaries  using  a 
Baratron  pressure  transducer  to  monitor  the  upstream  pressure. 

H atoms  were  monitored  using  the  H + NO  glow  compared  to  the 
0 + NO  glow.  The  stable  gas  concentration  was  monitored  by 
observing  the  partial  pressure  of  the  gas  in  the  absence  of 

[ 

oxygen  atoms  (or  any  reactive  gas)  . 

The  radiation  observed  is  analyzed  in  several  ways.  First, 
a synthetic  spectrum  is  prepared  by  a computer  model  for  com- 
parison with  experimental  spectra  (see  Part  II) . 

* 

Second,  a computer  model  of  the  reaction  kinetics  (e.g.,  Krie- 

■ 

ger  and  Kummler,  1973)  is  used  to  analyze  the  dominant  processes 
in  the  contact  time  and  concentration  regimes  of  experimental 
interest.  Third,  the  computer  model  is  employed  to  predict  re- 
gimes in  which  elementary  reactions  can  be  isolated  and  studied 
in  isolation  or  in  small  groups. 


[ 


T ? 
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RESULTS : SPECIFIC  SYSTEMS 


A.  H + N02 

The  reaction  H + NO^  -*■  OH^  + NO  was  chosen  as  a prototype 
reaction  to  study  because  it  is  the  classic  OH  generation  mech- 
anism for  flow  tube  work  and  therefore  must  be  quantified  be- 
fore other  OH  reactions  can  be  studied.  It  was  previously 
(Kummler,  1969)  shown  that  OH^  is  the  likely  product  of  that 
reaction  on  the  basis  of  kinetics.  However,  the  quantum  yield 

of  OH^  as  determined  from  absorption  measurements  by  Del  Greco 

and  Kaufman  (1962)  was  found  to  be  low.  The  study  conducted 
in  this  report  was  aimed  at  re-examining  the  quantum  yield 
via  infrared  emission  studies. 

The  procedure  followed  in  this  experiment  consisted  of 
calibrating  an  ambient  PbS  detector,  chopper,  filter  and 
PAR/HR8  lock-in  amplification  system  using  the  infrared  data 
of  Stair  and  Kennealy  (1967),  Kaufman  (1974),  and  Vanpee  (1971), 
for  the  photon  yield  of  the  reaction 

NO  + 0 - N02  + hv  . 

This  reaction  provides  an  actinometric  standard  at  2.7t  for 
our  system  once  the  NO  concentration  is  measured  via  an  MKS 
Baratron  and  the  0 atom  concentration  is  found  by  N02  titration. 
The  H atom  concentration  is  estimated  using  the  reaction  of 

H + NO  HNO*  -*•  HNO  + hv  . 
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The  observed  spectrum  from  H + NO  is  illustrated  in  Figure  3 
for  the  visible  region  of  the  spectrum. 

The  low  resolution  infrared  region  of  the  spectrum  for 
H + N02  is  illustrated  in  Figure  4 and  has  been  discussed  in 
previous  work  by  Kummler,  et  al,  (1973).  In  this  latter  work, 
a flow  velocity  of  80  m/sec  is  employed  in  order  to  minimize 
the  contact  time  (about  1 m sec)  following  the  introduction 
of  NC>2  to  the  discharged  stream.  Since  the  rate  constant 
for  H + N02  at  room  temperature  is  4.9  x 10-11  cc/sec,  and 
the  rate  constant  for  OH  + OH  + H20  + 0 is  2.5  x 10-12  cc/sec, 
we  have  suggested  that  the  emission  at  this  short  contact 
time  with  low  NC>2  concentrations  (4  - 5y)  is  primarily  OH.^ 
Thus,  in  Figure  4 we  see  the  0-1  and  1-2  transitions  near  2.7y 
and  the  0-2  transition  near  1.5y.  These  are  the  only  transi- 
tions possible  because  of  the  energetics  of  the  H + N02 
reaction  (£H2g=  -2843  kcal,  Av  < 2 for  OH).  Using  the  NO  + 0 
glow  to  calibrate  the  relative  intensities  of  the  1.5  to  2.7y 
radiation,  we  obtain  an  intensity  ratio  of  0.445  compared  to 
0 . 44 ±0 . 03  obtained  by  Murphy  (1971).  Without  knowledge  of 
the  detailed  OH^  destruction  reaction  rates  in  this  reaction 
system,  we  cannot  precisely  determine  the  quantum  yield  of  OH^, 
but  based  >n  ground  state  chemistry  we  have  concluded  that 
Kaufman's  estimate  appears  reasonable,  i.e.  the  quantum 
yield  is  much  less  than  0.10.  The  purpose  of  this  work 
is  to  re-examine  this  question  using  higher  resolution  spectro- 
scopy and  a detailed  examination  of  the  system  kinetics. 
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Our  initial  hypotheses  for  the  H + N02  system  were  formu- 
lated on  the  basis  of  low  resolution  spectra  illustrated  in 
Figure  5.  Based  on  the  work  of  Kadaj  (1973' , it  appeared  that 
O+C^H^  gave  vibrationally  excited  OH  which  radiated  at  2.7y 
via 

CHO  + 0 -*•  CO  + OH* 

and  at  1.5y  in  the  overtone.  The  reaction  of  CHO  + 0 is  exo- 
thermic by  79  kcal  and  can  therefore  produce  hot  bands  in  both 
the  fundamental  as  well  as  the  overtone.  Hot  band  detection 
in  the  fundamental  cannot  be  carried  out  with  great  sensitivity 
because  the  PbS  detector  loses  sensitivity  beyond  3y.  The 
reaction  of  H + NO2  was  presumed  to  yield  OH  fundamental  and 
overtone  radiation  in  the  primary  step,  although  the  hot  bands 
are  not  probable  since  this  reaction  is  only  exothermic  by 
2h . 4 kcal,  resulting  in  a maximum  excitation  of  v = 2 in  the 
product  OH.  No  bands  at  2.3y  such  as  those  present  in  the 
0 + C2H^  spectrum  are  found  which  is  consistent  with  the 
conclusion  that  the  2.3y  band  from  0 + C-H,  was  CO  overtone 
radiation.  As  shown  in  Figure  5,  the  0 + N2H^  spectrum  appears 
to  be  significantly  different  and  on  the  basis  of  the  expected 
kinetics,  it  is  probable  that  the  2.7y  radiation  is  due  to  H20 
( \>1  +V3)  . The  presence  of  a 1.9y  band  is  also  consistent  with 
the  presence  of  water  emission  in  the  0 + N2H^  reaction  (vida  infra) . 
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FIGURE  5 
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A comparison  of  low  resolution  spectra  originatir. 
from  0 + , 0+C2H4,  and  H + NOj  under  comparabi 

conditions . 
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With  these  previous  results  in  mind,  we  have  re-examined 
the  H+N02  system.  In  Figure  6,  we  have  illustrated  a higher 
resolution  (X/AX  = 300)  spectrum  resulting  from  that  chemilu- 
minescent reaction.  Similar  resolution  spectra  have  been  ta- 
ken of  0+C^.H,  and  no  structure  such  as  that  found  in  Figure 
6 has  appeared.  We  must  therefore  conclude  that  the  emission 
in  the  0+C2H4  system  is  not,  as  we  had  previously  concluded, 
emitted  from  the  same  molecule  albeit  with  higher  excitation, 
as  that  occurring  in  the  H+N02  system.  Murphy  (1974)  has 
taken  high  resolution  interferometric  observations  of  the 
H + NO 2 system  and  concluded  that  both  NO  and  OH,  but  not  H20 
emission  is  present  in  the  chemiluminescent  spectrum.  Our 
observations  are  consistent  with  OH^  being  the  major  emission 
source  in  the  2.7y  region  for  H + NC>2.  From  the  structure 
observed  in  this  OH^  spectra,  we  believe  that  OH^  is  not 

the  emitting  species  in  the  0 + C2H4  system  at  2.7y. 

To  assist  in  drawing  these  conclusions  from  our  spectra, 
synthetic  spectra  of  NO  and  OH  were  drawn  via  the  Calcomp 
plotter  (see  Part  II).  These  are  shown  in  Figures  7 and  8. 
Precise  agreement  with  the  experimental  spectra  is  lacking  for 
several  possible  reasons.  First,  there  is  an  uncertainty  of 
AX  = 0.002y  in  the  wavelength  position  of  the  experimental 
spectra  due  to  gearing  of  the  monochrometer  and  recorder  drive 
trains.  Second,  the  chemiluminescent  reaction  of  H + N02 
does  not  necessarily  produce  a Boltzmann  distribution  of 
vibrational  levels;  an  assumption  inherent  in  the  models 
illustrated  in  Figures  7 and  8. 
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Figure  8:  Artificial  Fundamental  Spectrum  of  OH  using 

Spectroscopic  Parameters  from  Herzberg. 


For  this  latter  reason,  we  obtained  the  high  resolution  OF 
spectrum  from  the  H + 0^  reaction  from  Murphy  (1971) . Since 
our  measurements  were  made  with  a relatively  low  resolution 
monochromator  (i.e.  1000y)  , we  integrated  Murphy's  line  by 
line  measurements  over  a 1000W  slit  formation  to  give  the  low 
resolution  spectrum  shown  in  Figure  9.  The  prominent  Q-branch 
transitions  of  Murphy's  integrated  measurements  agree  well 
with  the  prominent  features  shown  in  Figure  6.  Due  to  the 
low  exothermicity  of  the  H + N02  reaction  as  compared  to  the 
H + 03  reaction,  the  peak  distribution  is  different  between 
the  two  reactions.  However,  the  agreement  between  the  Q-branch 
features  leads  us  to  conclude  that  OH^  is  the  major  radiator 
at  2 . 7y  from  H + NC>2.  The  contribution  of  the  NO  overtone 
to  the  2 . 7p  emission  is,  at  present,  unresolved  since  only 
preliminary  measurements  have  been  made  in  the  5-6. 5y  region 
and  these  are  complicated  by  the  possible  presence  of  II20. 


B.  0 + N2H4 


\ 


N u . H~0  was  heated  under  vacuo  and  introduced  downstream 

b *1  b 

of  the  0 discharge.  A low  resolution  profile  of  the  resulting 
0 + N2H4  was  9^ven  F^9ure  5.  A higher  resolution  spectrum 
is  given  in  Figure  10.  Where  the  NO  overtone  emission  predomi- 
nates, according  to  the  NO  synthetic  spectrum  illustrated  in 
Figure  7,  no  emission  is  found  in  the  0 + N2H4  spectrum.  On 
the  other  hand,  a new  band  between  2.5  and  2.6y  appears  in 
0 + N2H4>  This  peak  caused  the  broader  appearance  of  the  0 + N2H4 
band  shape  relative  to  0 + C2H4  and  H + N02  in  the  comparison 
of  Figure  5.  However,  the  prominent  peaks  at  2.76  and  2.9y 
and  the  regularly  spaced  peaks  between  2.9  and  3.2y  do  closely 
resemble  H + N02  although  the  intense  :y  fall  off  with  increasing 
wavelength  is  much  less  than  for  H + N02-  In  fact,  the  intensity 
is  undoubtedly  increasing  near  3.2y  because  the  spectral  response 
of  the  PbS  tube  and  the  efficiency  of  the  monochrometer  (2.1y 
grating  blaze)  are  both  falling  off  rapidly  in  that  spectral 

region. 

The  reaction 

0 + N2H4  -*•  N2H2  + H20 

is  89.24  kcal  exothermic.  The  reaction 

0 + N2H4  -*■  N2H3  + OH 

is  only  36.12  kcal  exothermic  and  is  reported  by  Foner  and  Hudson 
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(1970)  to  be  only  l/25th  as  fast  as  the  formation  reaction. 

Thus,  the  branch  producing  OH  has  a rate  constant  of  about 
-12 

10  cc/sec  and  is  only  as  exothermic  as  H + NC>2  which  is  twenty 

times  as  fast.  Thus,  it  is  difficult  to  explain  the  apparent 
similarity  in  wavelength  and  band  spacing  of  0+N2H4  and  H + NC>2 
between  2.9  and  3.2y.  In  this  case  the  intensity  changes  appear 
to  be  more  significant  than  the  wavelength  positions.  The  fact 
that  emission  is  observed  out  to  2.5y  strongly  suggests  that 
H20  is  a major  contributor  to  the  2.7y  emission  as  shown  in 
Figure  10.  Kinetic  modeling  of  this  reaction  system  has  not 
yet  been  completed  to  further  verify  these  conclusions. 
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C.  0 + NR. 


In  order  to  evaluate  the  effect  of  substitution  on  amines, 
oxidation  by  0 atoms,  a series  of  experiments  using  tertiary 
amines  were  run.  Trimethyl,  triethyl,  and  tripropyl  amines  were 
introduced  by  beating  the  liquids  at  reduced  pressure  and  bubbling 
argon  through  the  heated  liquid.  Because  of  the  greatly  varying 
vapor  pressures  of  these  compounds  it  was  difficult  to  control 
the  partial  pressures  to  preset  values.  Nonetheless,  from  the 
observed  spectra  several  conclusions  can  be  drawn.  In  the  near 
ultraviolet  and  visible  spectral  region,  as  shown  in  Figure  11, 
the  most  prominent  peak  is  the  OH (A  - X)  electronic  emission 
at  3064  A (and  the  second  order  at  6128  A) . Also  observed  are 
the  CH (A  -*•  X)  at  4 312  A and  an  enhancement  of  the  7600  A 
region  which  also  appears  in  the  background  from  the  discharged 
0^.  Comparing  this  to  our  0 + NH^  spectrum  the  likely  radiator 
is  HNO*.  It  is  evident  from  a comparison  of  these  tertiary 
amine  + 0 spectra  with  0 + HC  spectra  taken  from  previous  work 
(Krieger  and  Kummler,  1973)  in  Figure  12  ■ .at  the  hydrocarbon 
fragments  play  a dominant  role  in  the  chemiluminescent  process. 

In  the  infrared  spectra  illustrated  in  Figure  13  we  also 
find  that  our  spectra  for  tertiary  amines  resemble  those  for 
hydrocarbon  + 0 illustrated  in  Figure  5.  This  fact  is  parti- 
cularly apparent  by  reference  to  the  2.5  - 2.6y  region  where 
the  substituted  amines  show  little  emission  in  contrast  tc 
hydrazine  but  in  agreement  with  0 + No  infrared  signature 

for  0 + NH^  is  observed  in  the  2.7y  region. 
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Figure  11:  Comparative  visible  chemiluminescence  from  the 

attack  of  0 atoms  on  several  substituted 
amines.  The  prominant  feature  at  3100  A is 
OH*,  the  second  order  of  which  is  at  6200  A. 
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L3:  Comparative  infrared  spectrum  for  the  0 atom 

attack  on  several  substituted  amines. 
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Introduction 


In  order  to  interpret  the  extensive  infrared 
measurements  described  in  the  previous  section,  an 
analytical  model  was  developed  to  predict  the  non- 
equilibrium band  spectra  associated  with  several 
diatomic  molecules  of  interest.  This  section  describes 
this  analytical  model  and  the  "'omputv.r  code  used  to 
generate  artificial  infrared  spectra*.  The  code  uses 
anharmonic  energy  parameters  and  line  strengths  as 
obtained  from  Herzberg^^.  The  vibratjonal  distributions 
for  the  diatomic  species  can  be  assumed  as  either 
Roltzmann  or  input  in  any  form  desired.  The  rotational 
distributions  are  assumed  to  be  Roltzmann  which  is 
consistent  with  the  experiments  above. 

The  program  performs  a line-by-line  calculation 
for  a predetermined  number  of  hot  bands,  then  integrates 
the  compiled  intensity  versus  wavelength  values  over  an 
arbitrary  slit  function.  Figure  1 show?  a block  diagram 
of  the  program.  Input  data  to  the  MAIN  routine  includes 
the  number  of  vibrational  levels  to  be  cons idered  ,N ; the 
maximum  vibrational  level  included , IVMAX ; the  initial 
vibrational  level, IV;  the  v transition  to  be  considered, 
IVP;  a Q-branch  key, LAMBDA;  the  vibrational  and  rotational 
constants , WEI , XE1 ,YE1 , BE  and  AE;  an  appropriate  set  of 

*preliminary  code  written  by  R.  Marriott. 
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1 Assigns  constant  values 
? Vibrational  level-populations 

3 Vibrational  energies 

4 Calls  SPCTRM  I 


1 Rotational  level  populations' 

2 Total  line  energies  I 

3 Intensity  J. 

4 Calls  AMLGMT  J 


1 Orders  data  I 

2 Integrates  with  slit  function 

3 Calls  SLIT  ^ 


1 Calculates  slit  function 


PLOT 


1 Plots  graph  intensity  vs 
wavelength 


Figure  1:  Bl^ck  Diagram  of  the  Spectra  Code 
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Einstein  A-Coe£f icients ,A(I) ; the  rotational  and 
vibrational  temperatures ,TE  and  TV;  the  slit  width, 

S;  and  a parameter ,DE , used  in  controlling  the  movement  of 
the  slit  function  across  the  line-by-line  spectra. 

The  MAIN  routine  calculates  the  vibrational  energies  and 
population  distributions  then  calls  subroutine  SPCTRM. 

j 

The  maximum  J value  required  for  the  calculation  is 

1 

calculated  in  SPCTRM  followed  by  the  line  energies  for 

the  P , R , and  Q branches.  The  emission  intensities  for  I 

each  of  the  permitted  transitions  are  then  calculated 

using  the  appropriate  A-coeff icients . Subroutine  » 

J 

AMLGMT  is  then  called  to  order  the  line  energies  in 
order  of  decreasing  wavelength.  The  corresponding  line 
intensities  are  then  integrated  over  a triangular  or 
Gaussian  slit  function  in  energy  increments  of  DE.  The 
spectral  resolution  of  the  calculation  is  governed  by 
the  input  slit  function  width  at  half  heighth,S.  Subroutine 
PLOT  then  takes  the  assembled  intensity  versus  wavelength 
data  and  generates  a computer  plot  of  the  artificial 
spectra . 

A compilation  of  non-equilibrium  spectra  generated  by 
the  code  is  given  in  the  following  Figures.  Verification  of 
the  computer  generated  spectra  for  the  case  of  CO  has 
been  done  and  shows  very  good  agreement  with  experimental 
measurements . 
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Comparison  of  the  Artificial  Spectra 
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with  Experimental  Infrared  Spectra 
of  CO  Overtone. 
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Effect  of  the  DE  value  of  the  Artificial 
Spectra  of  CO  for  stated  conditions. 
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appendix  a 


A copy  of  the  program  that  generates  the  infrared 
data  is  contained  in  Appendix  A. 


DIMENSION 
DIM  EN  S ION 
DIMFNSION 


PN  (5?)  ,TF(50,5P)  , PE  (50,  50)  , FEL  <500C).F3L  (5000) 
FELP  (5C0C)  ,FSLP(5PrC) 

A ( 20)  ,EN  (1C)  , Til  (20)  ,TI2  (20)  ,SLI  (9) 


KF =0 

RFAD(  5,  11)  (Til  (I)  ,1=1,4) 

READ  (5,11)  (T 12  (I)  ,1*1,4) 

F EA  D(  5, 1 1)  ( SLI  ( I)  ,1  = 1,2) 

11  FORMAT (6 AU) 

5 READ (5,  6,  END  = 12 5)  N,IVMAX,IV,I  VP, LAMBDA 
WRI  TS  (8  ,fi)  N , IV  MAX,  IV  , IV  P,  LAMBDA 

6 FORM  AT  ( 1 014) 

READ  (5,13)  WF1  ,XE1  ,YE1,  BE,  AE 
HR  ITF.  ( 8 , 1 3)  WEI  ,XE1  , YE1  , BE  , AE 
READ  (5,13)  (A  (I)  ,1=1,  N) 

WRITE  (8,13)  (A  (I)  ,1  = 1,  N) 

READ(5,13)  TF, ,TV,  DE,S 
WR ITF  (9,13)  TE  , TV  ,DE  , S 
13  FOPMAT(8Fln.  U) 

CM=1. 26E+1 4 
ILOOP=I VM AX+IV  P 

ZER  1=0. 5* HE  1-0. 25*XE1*C.125*YF1 
DO  20  J=1  ,1C 
DO  20  K=1,  10 
TR ( J , K)  =0  • 

20  RE(J,K)=D. 

EN  (1 ) =0 

DO  4C  J =1 , I LOOP 
R=J-.5 

40  EN (J)  =WF 1 *R“ XF 1 *P **2  + YEl *R**3-ZER1 
DO  5*  J=1 , IVMAX 


50 

60 

70 

80 

80 


15 

80 


JJ=J+IVP 

RE(J,JJ)  = EK(JJ)-EH(J) 

0 = 0. 

DO  60  J=1 , 1 LOOP 

Q = F.XP  (-EN  ( J)  *1 . 43886/TV) +0 

SUM=-I LOOP 

DO  70  J= 1 , ILOOP 

SUM=  SOM*  J 

DO  80  J=1, ILOOP 

PN  (J)  =CM/I  LOOP 

PN (J)  =CM/SUM* J-CM/SUM 

PN(J)  =-CM/SUM*J*CM/SUM*ILOOP 

PH  (J ) =CM*  EXP  (-EN  (J)  *1.  4 3886/TV)  /Q 

PN(  1)  =CM/0 

WRITE(10  ,15)  J,PN  (J) 

FORM  AT  ( II , E1C  . 3) 


CONTINUE 

DATA  A /0  . , C . , 0 . , 

DATA  A/0.  ,0.  ,0.  , 


0.,1  . ,4.55,12.4,28.4  ,50.  ,104./ 

1.0,3.77, 8. 9 4, .16-^,  28.  3,  4 2.7,  57.8/ 


no  o o o n 


DATA  A/0., 
DATA  A/n • , 
DATA  A/C. , 
DATA  A/C., 


''..1.0'*,  2. 8,5.2  .8.  1,1 1.1,14.1,  16.  6,  18.6/ 

1.C,  1.  36,1 .26, .878, .3 9,. 0474  ,.0753/ 

C.  ,1.0,  3. 0,5. 9 6, 9. 89,  14.  74,20.48,  27.03,34.16/ 
1.0,1.93,2.78,3.53,4.21 .4.81 ,5.35,5.74 ,6.23/ 


R AD  = 1 .r 

DO  100  J = 1 , IVMA  X 
TR  (J,  J*IVP)  = A (J  + IVP)  *PAD 
00  CONTINUE 

CALL  SPCTRM  (TE, LAMBDA,  IV MAX,  IV,  IV P,  DE,  PN,  TR.  S iJ,  K F , FSL, FEL, BE , AE , S) 
GO  TO  5 
125  CONTINUE 

WP  ITE  ( 6,  1)  KF 
1 FORMAT  (*KF= • ,14) 

WRITE  (6,150)  (FEL  (I)  ,FSL{I)  ,1=1  ,KF) 

150  FORMAT (/, 1P2E10.3) 

DO  ID  1=1, KF 


I I=K  F ♦ 1- 1 


FELP  (I)  =1.  E+04/FEL  (II) 

10  FSLP  ( I)  =PSL  (II) 

WRITE  (7,150)  (FEL  P (I)  , FSLP  (I)  , 1=  1,KF) 

CALL  GR  AP  H ( FFLP , FSLP ,KF) 

CALL  PLOT  (FELP,  FSLP,  KF,TI1,TI2,  SLI) 

STOP 

END 

SUBROUTINE  SPCTRM  (TE  , LA  MBD  A , I VM  A X , I V ,1 V P - DE  , PN  ,T  R,  RE,  KF,  FS  L,  FFL,  BE,  AE, 
DIMENSION  S L (5000  ) ,EL(50C0)  ,WSL  (500),  W EL  (500)  ,PN  (50)  , NT  (100) 

DIMENSION  PSL  (10,5000)  , PEL  ( 1 0 , 5000)  ,TR(5C,5C)  ,RE(50,50) 

DIMENSION  FSL(5C00) ,FEL(5000) 

AO=1 . 44/TE 
DO  1 Cn  I=IV , IV  MAX 
DO  107  IS=1,10C0 
SL  (IS)  *0.  ? 

107  EL  ( IS)  =C.d 
KSNT=9 


IVP1  = IVP*I 
BV=BE-AE* (1-0.5) 

BVP  = BE-AF*  (IVP1-C.  5) 

JD=  - 3 
A= AO*BV 
FJ  1 =EXP  (- A) 

EJ2=EJ1*EJ1 
AOR  = 1.0 
QP=1 .0 

JC=0  « 


JMAX='> 

101  *1MAX=JMAX+1 
JL=JC*1 
JC2=JC  + 1 


DO  1 C 2 IJ=1 ,JC 
102  A0R=A0P*EJ1 
AQR1=  AOR* JC2 

5 pO  R M A T ( 1 AQR1  = ' , 1PE10.3,'  JMAX='  ,15) 

QP=QR+ AQR1 

IF ( AQR1/0P-  1.E-  3)  10  3,103,101 
1C3  EBL=RE  (I  ,IVP1 ) 

C3  = PN  ( T V P 1 ) *TR  (I,IVP1)  /OK*1 .98*35-23 
C WRITE  (6,200rj  JMAX  , FBL,  CB,  PN  ( IV  P 1)  ,TR  (I,  IVPl  ) 

2000  FORM  AT  ('  JMAX=»,I4,'  SBL=  • ,E  10. 3,  ' CB=',E10.3 
1,'  PN*',E10.3, • TR  = ' , E 1 0 . 3 ) 

DO  1 C 4 JP=1  ,2  H-25 
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J D=iT  1)  + 2 
EJ=1 . 0 
FFJ  = 1 . O/EJ  2 
WPITF  (6  , 1 2 C ) 

120  F0PHAT(//,‘  WSL ' , 8X, * WSL' ) 

DO  IC  S J= 1 , JHA  X 
EFJ=  EF  J*  EJ2 
JD1  =J  *0D 

WEL  (J)  = EBL+BVP*J*  (J-1 ) -BV*(JD1-  1)  *JD1 
EJ  = EJ  * ^ F J 
JDE  H= J+  JD1 - 1 
IF(JDFH)  11,12,11 

12  SJ=C 

GO  TO  13 

11  SJ=  (JDEIHLAHBDA)  * ( J DEM-L  AM  BD  A) 

SJ  = SJ /IDEM 

13  WSL  (0)  =CB*WEL  (J)  *S  J*EJ 

WR ITF ( 6 , 20  A c)  WEL  (J)  , WSL  ( J) 

2 COS  FOPHAT  (1P2E10.  3) 

10  S COH  TIM  HE 

CALL  AHLGHT  (O.C  , JMAX,WSL,WEL,SL,  EL,  K ENT,  S ) 

: WRITE  (7, 20  1C)  KENT 

20 10  FORK AT ( ' KENT  = ' , 14) 

104  CONTINUE 

IF(LAHBDA)  113,112,113 

113  EJ=1.n 

WP  ITS  ( 6,  1 2 C) 

EF  J=1 .C/EJ2 
DO  114  .7  = 1, JM AX 
EFJ=EFJ*  EJ2 

WEL  (J)  =FBL*BVP*J*  (0-1)-BV*J*  ( J-1 ) ♦ (B  V-BV  P)  *LAM3DA  * * X 
EJ=EJ*EFJ 

SJ=  (2*  J + 1)  *LAHBDA*LAMBDA 
SJ=SJ/(J*(  J+1)  ) 

WSL  (J)  =CB*  WEL  (0)  *S J*EJ 
WRI7E(  6, 200  5)  WEL  (J)  ,WSL(J) 

114  CONTINUE 

CALL  AMLGHT(0.0,JMA.  X,WSL,WEL,SL,EL,KENT,S) 

: WRITE (7 ,20  IF)  KENT 

112  KT  (I)  =KENT 

DO  1 C 6 KS  = 1.,  KENT 
FSL  (I  ,KS)  = SL  (KS) 

FEL  ( I,  K S)  =F  L (KS) 

106  CONTINUE 
WRTTE{6,  1) 

1 FOPHAT (//,'  PEL* ,8X, ' RSL' ) 

WRITE  (6,  2015)  (REL  (I  ,1  J)  ,RSL(I,IJ)  ,IJ=1,KENT) 

2015  FOPKA?( 1P2E10.3) 

CALL  AKLGHT  (DE,  K ENT  , S L , EL,  FSL,  FEL,  KF,  S) 

: WP  ITF  ( 6 , 2) 

2 FOPHAT  (//,'  FEL'  ,8X,  ' FSL'  ) 

WPITF  (6, 2C  15)  (FEL  (J)  ,FSI,(J)  ,J=1,KF) 

: WRITE  (7 , 20  2C)  KF 

2 02  0 FOPHAT  {'  KF=',I4) 

100  CONTINUE 
RETUPN 
FND 

SUBROUTINE  AHLGHT  (D,  J , WS , W E,  S , E, K,  SL) 

DIHEN  SION  WS  (50  00)  ,WF.  (5000)  , TABE  (5000)  ,TABS(500C)  ,S  (50^0) 
J1=K  11-26 


, E (5000) 


*5 


VNMlhM* 


r 

IF  (WE(J)-WE(I))  100,101,101 

1 

101 

NW  = 1 

2 

NM=  J 

.3 

N D=  1 

4 

GO  TO  102 

8 

100 

NW=J 

A 

NM=1 

7 

ND=-  1 

8 

102 

IF  (D)  1r 3, 103, 104 

q 

103 

IF  ( K)  1 25,  125,  126 

V 

125 

DO  127  1=1, J 

1 . 

S ( I)  =WS  (NW) 

0 

tm 

E (I ) = WE  (NW) 

3 

127 

NW=NW+ND 

u • • 

K=J 

8 

GO  TO  1000 

6 

126 

JCOUNT= J+K 

7 “ 

IC  = 0 

8 

NC=0 

9 

K C=JCO  UN  T 

1 42 

ETEST=  WE  (1  ) 

1 

IE=  1 

2 

NC=NC+1 

3 

DO  131  1=1, J 

4 

IF  (ETEST-WS  (I) ) 132,131,131 

5 

132 

IE=  I 

6 

ETEST=  WE (I) 

7 

131 

CONTINUE 

8 

DO  133  1=1  , K 

9 

IF ( ETEST- E ( I)  ) 144,133,133 

0 

1 44 

I E = -I 

1 

ETEST=E  ( I) 

2 

133 

CONTINUE 

3 

IF ( I J)  105,105,106 

4 

105 

IE=“IE 

5 

TABS(JCOUNT)  =S  (IE) 

A 

TA  BE  (JCOUNT)  = E (IE) 

7 

E ( I E)  =-  E ( IE) 

8 

GO  TO  107 

9 

106 

TABS  ( JCOUN  T)  =WS  (IE) 

c 

TABF (JCOIJNT)  = WE  (IE) 

1 

WE(  IE)  =-  WE  ( IE) 

2 

107 

IF  (JCOUNT-KC)  108,109,109 

3. 

103 

IF(  TAB E ( JCO  UN  T)  - TAP  E ( JC  OUNT+ 1 ) ) 109,141 

4 

141 

JC1 = JCOUNT* 1 

5 

TABS  ( JC  1)  = TABS(JCOUNT+1)  + TABS  (JCOUNT) 

ft 

IC=IC-1 

7 

GO  TO  110 

8 

1 09 

JCOUNT= JCOUNT -1 

19 

1 10 

IF(KC-NC)  142,  143,142 

'0 

143 

K=  KC+ IC 

! 1 

DO  14  5 I = 1 , K 

!2 

E (I ) = T A B E ( I- IC) 

145 

S(I)  =TABS(  I-  IC) 

’4 

DO  146  1 = 1 , J 

>5 

1 46 

WE(  I)  =-  WE  ( I) 

>6 

GO  TO  1000 

>7 

104 

L E = 2 

>8 

IF (K)  111,111,112 

>9 

111 

RN=WE(NK)  /D 

' TT 
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HI 


RN 
= 1 
F.O=N*D 
FS?ORF,=  EO 

NSTF.PS=  (WE  (NM)-WE(NW)  ) /D  + 1 
J 1=NSTEPS 
EMAX=FO+NSTEPS*D 
GO  TO  114 

IF( WE  (NW)-EO)  115, 113,113 

N =N  - 1 

EO=EO-n 

NSTEPS=N  STEP  S + 1 
ESTOPE-  EO 
GO  TO  112 

IF  (WE  (MM) -EMAX)  114,114,  116 

N STEP  S=N S^EPS  + I 

EMAX=SMAX+D 

GO  TO  113 

DO  117  IE=1  , NS7EPS 

TABE(  IF)  =E  STORE 

TABS  (I E) =9  . C 

TAS=TALSS(  IE) 

CALL  SLIT  (WE,WS,T  AS,  ESTOPE,  J,SL) 

TABS  (TE)  =TA S 

GO  TO  (1  19 ,120)  ,LE 

DO  13  5 JE  = 1,J1 

IF  (E  (OF.)  -ESTORE)  135,137,  135 

TA3S  ( IF)  = TA  B S ( IE)  +S(JE) 

GO  TO  119 
CONTINUE 

WRITE(6, 1)  ESTORE, IE 

1 FORMAT  (/'  ESTORF.=  ’ , E15.  8,  ' IE=*  , 14) 

WRITE  (6,2) 

2 F 0 R M A T ( 1 E ( JE)  ') 

WRITE  (6, 3)  (E(JF)  , J E= 1 , NST EPS) 

3 FORM  AT  ( IX,  El  5.8) 

ESTOPE=ESTOPE+D 
CONTINUE 
ESTOFE=  EO 

DO  121  IE  = 1 ,N  STEPS 
E (IE)  = TABE  (IE) 

S ( I E)  =TA  B S ( IE) 

K= NST EPS 
3 RETURN 
END 

SUBROUTINE  SLIT (TE , TS , TSUM ,E STORE ,NENT ,SL ) 
DIMENSION  TE(50PO)  ,TS  (5CC0) 

WR  ITE  ( 10,10)  ESTORE 
10  FORMAT  ( El  0 • 3) 

SLT  = SL/10CC.*13.2E-7 
DLAM=SLT*  ESTORE*  ESTORE 

DO  ICO  1=1, PENT 
SFU  N = 0 . 

A=A  BS ( ESTOP  E- TE  (I) ) 

IF  (A.  LE.  DLAM)  S FU  N=  1 - A/DL  AM 
X = ( ESTORE- TF(  I)  ) * (ESTOPE- TE  (I)  ) 

IF  ( ABS  ( X)  . GT.93 . 54)  GO  TO  9° 

SPUN  = EXP  (-X/3B0) 

TSUM=TSUM  + SFUN*TS  (I) 

GO  TO  ICC  t t _ oo 


' u.,.^1  II I Xi J^II  1, 


■ ■ 1 "W1 


» W ■U,^W  1 -WVJ- 


C 

c 


FILE 


99  SFIJN  = 1C  .f‘E-10 

TSU  M=T SUM  + SFUN  *TS  (I) 

100  CONTINUE 
RETURN 
END 

SUBROUTINE  PLOT  (FEI, , F SL  , KP  , T1 1 ,TI2,SLI) 

DIMENSION  FFL  (500C ) , FSL(5CCC)  , Til  (20)  ,TI2(20) ,SLI  (8) 

CALL  PLTSIZ(I.O) 

X M I N = 2 . 50 
DX  = C.  1 

CALL  PSCALE(R.0,1 . 0,  XMTN,DX,FEL  (1 ) ,..P,1 ) 

CALL  PSCALF  (8 . 0 , 1 . 0 , Y tl IN,  DY  , FSL  ( 1)  , K F,  1) 

CALL  PLTOFS(XKIN,DX,YMIN,DY,1 .0,1 .0) 

CALL  PAXIS  (1.0,1.  C,'  WAVELENGTH  • , - 1 3, 8 . 0,  0 . 0,  XM  IN  , DX  , 1 . 0) 
CALL  PAXIS(1.9,  1.0,  ’INTENSITY  ( A.  0.)  * , 1 7 , 8 . 0 ,90 . 0 , Y MIN,  DY , 1 . 0) 
CALL  PLINF(1'SL(1)  , FSL  (1)  , KF,  1,0,  0,1) 

CL=  PSYMI.N  (0.2,16) 

CALL  PSYMB ( 5. 0-CL/2.f, 10.0,0.2, Til  (1)  ,0.0  ,16) 

CALL  PS  Y MB  (5 , O - PS  Y MLN  (0 . 15,  16)/2.,9,5,  0.15,TI2(1)  ,0.0,16) 

CALL  PSYHB  (5.0-PSYHLN  (C  . 1 5,8) /2 . 0 ,9. 1 ,0.  15,SLI(1)  ,0.0,8) 

CALL  PLTEND 

RETURN 

END 
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Introduction 

This  report  presents  the  results  of  vibrational  relaxation  in  a core 
flow  field  of  a rocket  engine.  Of  particular  interest  is  the  extent  of 
vibrational  freezing  in  a set  of  combustion  products  undergoing  rapid 
expansion.  The  relaxation  model  used  has  been  presented  elsewhere^)  but 
generally  includes  vibration-translation  (VT)  and  vibration-vibration  (VV) 
energy  transfer  processes  coupling  together  Hg,  Ng,  CO,  OH,  five  modes  of 
HgO  (100,010,020,001,  and  000)  and  eight  modes  of  COg  (100,010,020,030, 
001,101,021,  and  000).  The  model,  termed  the  "Berkeley"  model  after  its  site 
of  development,  contains  reactions  including  0 and  H which  are  not  present 
in  this  flow  field  and  state  of  the  art  rate  coefficient  data. 

The  flow  field  used  in  this  calculation  is  modeled  after  the  core 
flow  field  of  a TITAN  II  rocket  engine.  The  eight  representative  stream- 
lines were  supplied  by  Boynton^2)  and  feature  an  LTE  assumption  at  the  nozzle 
exit  plane  and  include  the  effects  of  film  cooling.  This  latter  effect 
produces  fuel  product  enrichment  (Hg,  CO  and  COg)  in  the  outermost 
streamlines  as  a result  of  the  fuel  rich  flow  used  to  cool  the  surface  of 
the  exit  nozzle. 

The  calculation  presented  here  uncouples  the  energy  transfer  chemistry 
from  the  flow  field.  This  assumption  is  generally  justified  since  the  heat 
capacity  of  the  flow  gases  is  only  weakly  dependent  on  the  extent  of 
vibrational  freezing  of  the  diatomic  species  present  in  the  flow  (vida  infra) 
The  general  conclusions  evident  from  this  relaxation  calculation  are: 
a)  The  diatomic  species  tend  to  freeze  out  early  in  the  expansion; 
near  exit  conditions  (900°K)  for  the  outermost  streamlines  while 
somewhat  lower  for  the  more  dense  inner  streamlines. 


b)  In  the  outermost  streamline,  all  modes  tend  to  freeze  near 
the  nozzle  exit  condition, 

c)  The  bending  modes  of  HgO  and  COg  tend  to  remain  in  LTE  with  the 

flow  for  all  except  tne  outermost  streamline. 

d)  The  presence  of  vibrationally  frozen  Ng  causes  COg(OOI)  to  be  above  LTE 

conditions  due  to  the  rapid  VV  coupling  which  in  turn  causes 
COg(lOl)  to  be  above  LTE  conditions. 

e)  Vibrationally  excited  Hg  is  predicted  throughout  the  core  flow  field 

but  it  is  not  yst  clear  wha-':  effect  this  will  have  on  mixing  layer  chemistry. 

The  following  sections  present  tne  details  of  the  calculation  and 
the  substance  behind  the  conclusions. 

Flow  Field 

(2) 

The  core  flow  field  was  obtained  from  a calculation  of  Boynton'  ' usino 

the  MULTITUBE  program.  The  calculation  used  as  initial  conditions  nonuniform 

(2) 

exit  plane  conditions  calculated  by  Aerospace  which  encorporatec  both  film 
cooling  and  nozzle  contour  effects.  An  inviscld  vacuum  expansion  frcm  these 
nonuniform  exit  conditions  was  performed  on  the  MULTITUBE  code^The  specific 
temperature  and  pressure  histories  on  eight  streamlines  across  the  vacuum 
expansion  are  shown  in  Figures  1 and  2.  Several  features  should  be  noted 
in  these  streamlines.  The  centerline  streamline  shows  recompression  as  ?. 
result  of  the  nozzle  contour  followed  by  a rapid  expansion  region.  The 
outermost  streamline  shows  the  most  rapid  expansion,  but,  of  course, 
represents  the  smallest  flow  mass  fraction.  And  the  variation  in  the  exit 
plane  temperature  across  the  nozzle  plane  is  clearly  demonstrated.  This 
latter  effect  is  a result  of  combustion  heating  leading  to  increased 
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temperatures  and  film  coolant,  e.g.  uncombusted  flows,  near  the  nozzle  walls. 

Table  I lists  the  exit  plane  concentrations  used  in  these  relaxation 
calculations.  As  expected,  fuel  rich  species  are  present  in  higher  concen- 
trations in  the  outermost  streamlines  as  a result  of  the  film  coolants. 

This  variation  in  composition  across  the  exit  plane  enhances  the  possibility 
of  having  vibrational ly  excited  hydrogen  entering  the  mixing  layer  region 
of  the  flow. 


Relaxation  Model 

The  (nergy  transfer  model  used  in  these  relaxation  calculations 
include  vit ration-vibration  and  vibration-rotation-translation  processes 
for  the  major  vibrational  modes  of  H^,  Ng,  CO,  HgO  and  COg.  Table  II 
gives  the  vibrational  levels  included  in  the  model  together  with  the  energy 
spacing. 

The  rate  coefficients  used  in  these  calculations  have  been  derived  from 
several  sources.  Primarily,  experimental  values  are  used  where  possible  but 
the  large  number  of  potentially  important  energy  transfer  processes 
necessitate  that  estimates  be  made  for  those  processes  for  which  no  experi- 
mental values  are  available.  Table  III  details  the  energy  transfer  processes 
included  in  the  model  together  with  the  references  regarding  the  rate 
coefficient  value.  Further  details  of  the  energy  transfer  processes  and 
rate  coefficient  estimates  can  be  found  elsewhere.^ 

Several  features  of  the  relaxation  model  should  be  indicated.  First, 
the  number  of  vibrational  levels  included  in  the  polyatomics,  i.e.  HgO  and  COg, 
are  fairly  large.  This  is  primarily  a result  of  interest  in  vibrational 
levels  of  the  polyatomics  which  give  rise  to  2.7p  radiation,  i.e.  COg (021 ) 
Cn2(101)  and  H20(001).  Secondly,  only  one  vibrational  level  is  included  in 
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Table  I 

Exit  Plane  Mole  Fractions 


Species 

Streamlines  1-4 

Streamlines  5-8 

H2 

0.04 

0.22 

h2o 

0.47 

0.33 

CO 

0.03 

0.10 

N2 

0.36 

0.32 

NO 

0.009 

0 0001 

co2 

0.094 

0.13 

1116  i 
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TABLE  III 


(2)  RELAXATION  MODEL 


Rate  Constant 

Vibration-Translation  Processes  Reference  (cc/molecule-sec) 


1) 

H2(v=1)  + 

H -*■  H2  + H 

(a) 

3.2"11 

exp(-1400/T) 

2) 

H2(v=1)  + 

0 -►  H2  + 0 

(b) 

3.0"13 

/t" 

3) 

H2(v=1)  + 

H2  - H2  + H2 

(c,d) 

1.5”  7 

exp(-139/T1/3) 

4) 

H2(v=1)  + 

h20  - H2  + h20 

(b,  e) 

1.5”  7 

exp(-139/T1/3) 

5) 

N2(v=1)  + 

H N2  + H 

(y) 

1.9”6 

exp(-164/T1/3) 

6) 

N2(v=1)  + 

0 «*■  N2  + 0 

(f /5/h) 

1.2”13 

1 exp(-23/T1/3) 

7) 

N2(v=1)  + 

H2  - N2  + H2 

(i) 

3.9”8 

exp  (-164/T1/3) 

8) 

N2(v=1)  + 

H20  -»■  N2  + H20 

<b,e) 

3.9”8 

exp  (-164/T1/3) 

9) 

OH (v=l)  f 

H + OH  + H 

(b) 

3-11 

10) 

OH (v=l)  + 

0 + OH  + 0 

(b) 

3”11 

11) 

OH (v=l)  + 

H2  + OH  + H2 

(b) 

l”11 

12) 

OH(v=l)  + 

H20  + OH  + H20 

(b) 

3”12 

13) 

H2O(010)  +H  -►  h2o(ooo)  +H 

(b) 

l”11 

14) 

H2O(010)  +0  H20(000)  +0 

(b) 

l”11 

15) 

H2O(010)  + h2  H20{000)  +h2 

(b,j) 

i-u 

16) 

H2O(010)  +H20  H20(000)  +h2o 

(k) 

r11 

17) 

H2O(001)  +H  H20{000)  + H 

(b) 

i”11 

18) 

H2O(001)  + 0 -►  H20(000)  +0 

(b) 

r11 

19) 

H20(001)  +H2  H20(000)  +h2 

(j) 

2”13 

20) 

H20(001)  +H20  -*■  H20(000)  +h2o 

(1) 

i”12 

1 1 1-8 


21) 

H20(001) 

+ H H20(020)  + H 

22) 

H20(001) 

+ 0 ■+  H20(020)  + 0 

23) 

H20(001) 

+ H2  ■*»  H20(020)  + H2 

24) 

H20(001) 

+ H20  -*■  H20(020)  +H20 

25) 

CO (v=l) 

+ H ■*>  CO  (v=0)  + H 

26) 

CO (v=l) 

+ 0 •*  CO  (v=0)  + 0 

27) 

CO (v=l) 

+ H2  CO  (v=0)  + H2 

28) 

CO (v=l) 

+ H20  ■*>  CO  (v=0)  + H20 

29) 

co2(oio) 

+ H -*•  C02(000)  + H 

30) 

CO2(010) 

+ 0 + C02(000)  + 0 

31) 

CO2(010) 

+ H2  C02(000)  + H2 

32) 

o 

o 

to 

o 

o 

+ H20  ->  C02(000)  +H20 

33) 

C02(010) 

+ C02  -*■  C02(000)  +C02 

34) 

C02(001) 

+ H ->  C02(000)  + H 

35) 

C02(001) 

+ 0 -*•  C02(000)  + 0 

36) 

C02(001) 

+ H2  C02(000)  + H2 

37) 

C02(001) 

+ H20  ■*  C02(000)  +II20 

38) 

C02(001) 

+ H -*■  C02(030)  + H 

39) 

C02(001) 

+ 0 -►  C02(030)  + 0 

40) 

C02(001) 

+ H2  C02(030)  + H2 

(b) 

r11 

(b) 

r11 

(b) 

r11 

(b) 

r11 

(O) 

5.3”7  exp (-119/T1/3) 

(P) 

l“8  exp(-96/T1/3) 

(q) 

9"9  exp (-119/T 1/3) 

(b) 

9”9  exp (-119/T 1/3) 

(b) 

3.8"12  exp (— 7/T1/3) 

(P) 

2.3"9  exp  (-77/T1/3) 

(r) 

7 . 6~12  exp  (-7/T1/3) 

(s) 

l"11 

(r) 

4.6"10  exp (-77/T1/3) 

(t) 

6 . 7~8  exp(-208/T1/3) 

(t) 

l"8  exp(-96/T1/3) 

(t) 

9-9  expt-lig/T173) 

(t) 

9-9  exp(-119/T1/3) 

(b) 

2.1317  T"6-34  exp (-3013/T 

-378.7/T1/3  + 1416/T2/3) 

(u) 

4 . 624  t“5,89  exp(-4223/T 

-672 . 7/T1/3  + 2683/T2/3) 

(v) 

4.2717t“6,34  exp (-3013/T 

- 378.7/T1/3+ 1416/T2/3) 

1 1 1-9 
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41) 

ce2(ooi)  +h2o  -►  co2(030)  +h2o 

(V) 

4.7l”40T4*54  exp(2081/T 

+ 454/T173-  1729/T^3) 

42) 

CO2(001)  +C02  C02(030)  + CO 2 

(V) 

9.1623t"5,89  exp (-4223/T 

- 672.7/T1/3-h  2683/T2/3) 

43) 

CO2(001)  + CO  •*>  CO 2 (030)  + CO 

(V) 

6.8723t“5,89  exp (-4223/T 

-672.7/T1/3+  2683/T2/3) 

44) 

CO2(001)  +N2  -»•  C02(030)  tN2 

(V) 

6.8723t”5*89  exp (-4223/T 

9 

-672.7/T1/3+  2683/T2/3) 

Vibration-Vibration  Processes 

45) 

H2(v»1)  + H20(000) 

H2(v*0)  + H2O(001) 

(m) 

2.5"14  T5/6  exp (-30.1/T173) 

46) 

H2(v»1)  + OH (v*0) 

H2(v=0)  + OH(v»l) 

(b) 

l"12 

47) 

N2(v«1)  + H2O(010) 

N2(v=0)  + H2O(001) 

(b) 

3~13 

48) 

N2(v=1)  + H20(000) 

*►  N2(v=0)  + H2O(010) 

(n) 

9.4”11  exp(-64/T1/3) 

49) 

H2O(010)  + H2Q(010) 

+ H20(020)  + H20(000) 

(b) 

l”12  /T 

50) 

CO2(001)  + H20(000) 

CO2(010)  + H2O(010) 

(b) 

r13 

III-10 

51) 


N2(v=1)  + C02(000) 

-*•  N2(v=0)  + C02(001)  (r) 

52)  N2(v=1)  + CO (v=0) 


Take  the  I 1.7l"6exp(-175.3/T1/3) 
1valult  |6.07"14exp(15.3/T1/3) 


-* 

■ N,(v=0)  + CO (v=l)  (v) 

Take 

the 

1.78  exp  (-210, 

largest 

-13  . 

value 

6.98  x exp (-25 

53) 

CO, (001) 

+ CO  (v=0) 

C02(000)  + CO(v=l) 

(b) 

1.56 

_11  exp  (-30 .1/T3 

54) 

CO2(101) 

+ H20(000) 

-*■  CO2(000)  + H2O(001) 

(b) 

2-13 

/t" 

55) 

CO2(021) 

+ H20(000) 

■*  C02(000)  + H2O(001) 

(b) 

2-13 

/t~ 

56) 

CO2(101) 

+ C02(000) 

-*•  CO2(100)  + CO2(001) 

(b) 

2-11 

/t" 

57) 

C02 (020) 

+ C02(000) 

-*•  CO2(010)  + CO2(010) 

(b) 

4'13 

/F 

,-11 

58) 

H20(020) 

+ H ■*>  H2O(010)  + H 

(w) 

2 

.-11 

59) 

H20(020) 

+ O H2O(010)  + 0 

(w) 

2 

.-11 

60) 

H20(020) 

+ H2  **  H2O(010)  + H2 

(w) 

2 

.-11 

61} 

H20(020) 

+ H20  H2O{010)  +h2o 

(w) 

2 

62) 

C02 (020) 

+ H -*■  CO2(010)  + H 

(w) 

7.6“ 

12  exp(-7/T1/3) 

63) 

C02 (020) 

+ O CO2(010)  + 0 

(w) 

4.6“ 

9 exp  [-77  Z^3) 

64) 

CO, (020) 

+ H2  -*■  CO2(010)  + H2 

(w) 

1.5“ 

’U  exp  (-7/T1^3 ) 

65) 

C02 (020) 

+ H20  CO2(010)  + H20 

(w) 

2-ll 

.1/3, 


III-ll 


66) 

o 

o 

to 

o 

to 

o 

+ co2  C02(010)  +co2 

67) 

CO ^ (030) 

+ 

H C02(020)  + H 

68) 

C02(030) 

+ 

0 CO 2 (020)  + 0 

69) 

C02(030) 

+ 

H2  C02(020)  + H2 

70) 

C02(030) 

+ 

H20  C02(020)  + H20 

71) 

C02(030) 

+ 

C02  + C02(020)  + C02 

72) 

H20(100)  and  H20(001)  assumed  in 
local  translational  temeprature. 

73) 

co2(ioo) 

+ 

M -►  C02(020)  + M 

(w) 

9.2-10 

exp(-77/T1/3) 

(w) 

1.2"11 

exp(-7/T 1//3 ) 

(w) 

6.9-9 

exp(-77/T1/3) 

(w) 

2.3'11 

exp(-7/T1/3 ) 

(w) 

3’11 

(w) 

1.4'9 

exp(-77/T1/3 ) 

(X) 
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(5)  RADIATIVE  LIFETIMES 


SPECIES 

OH(v»l) 

H20(100) 

h20(oio) 

H20(001) 

H20(020) 

C02(010) 

C02(001) 

CO2(101) 


t (sec) 

0.3 

0.14 

0.045 

0.014 

0.022 

0.33 

0.0025 

0.091 


the  diatomic  species  due  to  a general  lack  of  'rate  coefficient  data. 

As  well,  this  model  is  not  interested  in  anharmonic  vibration-vibration 
effects ^ which  further  justifies  the  inclusion  of  only  one  diatomic 
vibrational  level.  Due  to  the  Fermi  resonance  between  the  100  and  001 
levels  of  HgO,  the  relaxation  model  assumes  a quasi-equilibrium  between 
these  levels  at  the  local  thermodynamic  temperature.  Lastly,  radiation 
losses  are  included  in  the  model  although  for  the  streamlines  considered 
here  these  losses  are  small. 
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Results  and  Discussion 


The  results  of  relaxation  calculations  on  the  eight  streamliners 
are  shown  in  Figures  3 (centerline)  through  Figure  10(outermost  stream- 
line). The  vibrational  populations  are  given  in  terms  of  a vibrational 
temperature  defined  through 


i ■ 


• th 


where  and  N0  are  the  density  of  the  i level  and  lowest  level  density. 


• N 


respectively,  e..  is  the  energy  of  the  i—  level  and  Ty  is  the  vibrational 


temperature.  As  expected  f^om  the  temperature  and  pressure  profiles  of 
Figures  1 and  2,  the  outermost  streamline  shows  the  most  vibrational 
freezing.  The  expansion  is  sufficiently  rapid  along  this  streamline 
(Figure  10)  that  essentially  all  molecules  freeze  near  the  exit  conditions. 

For  streamlines  that  are  progressively  closer  to  centerline,  the  bending  modes 
of  CO^  and  H20  begin  to  equilibrate  with  the  rapidly  decaying  transitional 
temperature  but  diatomic  molecules  continue  to  freeze  at  vibrational 
temperatures  only  slightly  below  exit  temperatures. 

The  general  features  of  the  relaxation  process  are  evident  from  the 
Figures.  Except  in  the  outermost  streamlines  where  the  expansion  process 
is  very  rapid,  the  bending  modes  remain  near  local  thermodynamic 
equilibrium  (LTE).  The  stretching  modes  of  COg  and  HgO  generally  remain 
slightly  above  LTE  values  while  the  diatomic  modes  remain  frozen  at  the 
highest  values.  The  extent  of  vibrational  freezing  is  small  along  the 
centerline  streamline  where  the  density  and  temperature  are  relatively  high. 

The  important  energy  transfer  processes  which  control  the  relaxation 
processes  are  evident  from  a review  of  the  rates.  Table  Iv  presents 
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Species 

TE 
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CO 

OH* 

H20(100) 
H20(010) 
H20(020) 
H20(001) 
C02(100) 
C02(010) 
C02(020; 
CO3 (030) 
C02(001) 
C02(101) 
C02(021) 
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a representative  history  of  the  rates  and  rate  coefficients  for  stream- 
line 7 during  the  latter  phases  of  the  relaxation  process.  For  each 
of  the  diatomic  species,  CO,  H2  and  N2  the  major  loss  process  is  a 
VV  process  generally  involving  a polyatomic  species.  CO  tends  to  freeze 
at  '.he  highest  temperature  with  the  major  loss  process  being  the  VV 
pumping  of  C02  to  C02(001).  N2*  relaxes  somewhat  more  rapidly  than  CO 

due  primarily  to  the  more  rapid  coupling  of  N2  to  C02(001).  H2  freezes 
close  to  the  N2  value  with  the  major  loss  process  being  the  VV  pumping 
of  H20(001).  None  of  these  diatomic  molecules  have  Vf  loss  processes 
which  compete  with  the  VV  processes  in  magnitude  at  least  for  the  species 

mole  fractions  encountered  in  these  streamlines. 

Both  CO*  and  N?*  pump  C02(001)  so  this  polyatomic  mode  tends  to 
freeze  cut  substantially  above  the  bending  modes.  C02(101)  tends 
as  well  to  freeze  out  above  the  bending  modes  since  this  combination 
mode  is  formed  with  nearly  gas  kinetic  collision  figuring  from  C02(001) 
and  C02(100).  These  excited  C02  modes  give  rise  to  4.3y[C02(001 )] 
and  2.7y  [C02( 1 01 ) ] in  amounts  above  that  expected  from  LTE  values. 

Lastly,  these  calculations  suggest  that  H2  enters  the  mixing  layer 
region  of  the  plume  flow  field  as  a result  of  vibrational  freezing  in 
the  vacuum  core.  The  effect  this  H2*  has  on  chemiluminescent  processes 
in  the  mixing  layer  is  being  currently  investigated  particularly  through 

reactions  such  as 

0 + H2*  -*•  OH  + H . 

Theoretical  and  experimental  studies  on  the  role  of  vibrational  energy  in 
simple  combustion  reactions  such  as  this  will  clarify  the  importance  of 
vibrational  freezing  in  the  core  flow  field  on  plume  radiance. 
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Introduction 


Many  reactions  of  interest  in  the  upper  atmosphere  involve  collisions 
between  molecules  and  atomic  oxygen.  Theoretical  cross  sections  for 
vibrational  relaxation  processes  involving  H^O,  CO^,  CO,  Ng,  Hg  and  0 atoms 
have  been  determined  by  MarriCut.^  Many  of  these  theoretical  estimates 
are  the  only  available  results  to  use  in  model  non-equilibrium  calculations 
involving  these  molecules.  The  calculation  of  cross  sections  for  vibra- 
tional excitation  of  these  molecules  by  heavy  body  collisions  are  carried 
out  by  using  a numerical  close  coupling  solution  of  the  wave  equation  for 

the  colliding  system.  This  system  is  characterized  by  a breathing  sphere 

(21 

model  of  Schwartz,  Slawsky  and  Herzfeld'  ' to  represent  the  tarqet  molec- 
ular vibrator  and  by  an  empirical  spherically  symmetric  scattering  potential. 
Reliable  central  potentials  required  for  the  calculation  of  the  excita- 
tion cross  sections  were  not  always  available.  We  have,  therefore, 

calculated  these  potentials  from  scattering  cross  sections,  analogous  to 

(31 

potentials  determined  for  He-He  interaction.  ' Objections  to  the  approxi- 
mation of  anisotropic  interactions  by  spherically  symmetric  potentials  will 
be  considered  in  the  discussion  of  the  results.  Sinano^lu  has  suggested  an 

alternative  method  for  calculating  potential  parameters  in  the  additive 

(41 

force  approximation.  ' We  have  applied  this  model,  assuming  exponential 

interatomic  potentials,  to  the  collisions  of  0 atoms  with  homonuclear 

(M 

diatomics;  the  results  are  compared  with  those  of  Leonas.  w We  also  present 
the  raw  cross  section  data  and  scattering  configuration  so  that  other 
methods  can  be  applied  to  the  data  when  they  become  available. 


Experimental  Method 


Details  of  the  apparatus  have  been  given'  ' elsewhere  so  that  here 
we  discuss  only  the  scattering  geometry  and  problems  specific  to  this  wort. 

Ions  of  atomic  oxygen  are  produced  by  electron  impact  on  02»  pro- 
ceeding via  dissociative  attachment,  dissociative  ionization  or  ion  pair 
formation. ^ Either  positive  or  negative  ions  are  extracted  from  the 
source  and  momentum  selected  in  a transverse  magnetic  field.  0 ions  are 
unambiguously  distinguished  from  the  neighboring  OH  ions  by  calibration 
with  a Hall  probe.  Neutral  0 atoms  are  obtained  by  charge  transfer  with  02 
in  the  case  of  O',  and  with  Kr,  Xe,  N?  or  HgO  in  the  case  of  0+.  Belyaev, 
Kamyshov  and  Leonas  (BKL)  have  used  this  technique  to  measure  the  O-Ne 
potential  and  report  significantly  smaller  cross  sections  when  the  0 atoms 
are  produced  from  0+  and  Xe.^  Based  on  the  favorable  energy  defect, 
this  was  attributed  to  the  formation  of  metastable  0(  D)  in  the  charge 
transfer  reaction. 

The  scattering  configuration  has  cylindrical  symmetry  in  which  the 
entrance  hole  to  the  target  gas  cell  has  diam.  dg  = 0.0203  in.  while  the 
effective  diam.  of  the  Channel tron  detector  is  dQ  = 0.0362  in.  The  beam 
is  almost  parallel  giving  a rectangular  intensity  distribution  of  width  dg. 

The  exit  hole  of  the  scattering  cell  has  a diameter  of  0.025  in.  so  that  an 
optical  path  exists  between  the  detector  perimeter  and  all  points  in  the 
scattering  path.  The  distances  from  the  detector  to  the  entrance  and  exit 
holes  of  the  scattering  cell  are  respectively  1.431  in.  and  0.985  in.  Unless 
otherwise  stated,  the  experimental  incomplete  integral  cross  section  Hexp  shall 
be  understood  to  include  an  end  correction,  typically  2%,  to  take  account 
of  the  nonuniform  target  gas  density.  Thus  in  calculating  cross  sections 
for  comparison  with  experiment  we  may  assume  a uniform  density  inside  the 
cell,  and  zero  density  outside. 
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Results  and  analysis 

Qexp  was  measured  as  a function  of  lab.  energy  in  the  range  100-2500eV. 
The  interaction  potential  is  assumed  to  have  the  exponential  form  V(r)  = 

A exp(-ar).  By  introducing  reduced  variables  we  have  simplified  the  problem 
of  finding  A and  a to  a one-parameter  search  for  A,  a being  determined  by 
a least  squares  condition;  details  of  the  analysis  are  given  in  Appendix  A. 
The  results  are  given  in  Table  I.  With  the  exception  of  the  O-Ne  and  O-He 
potentials  (Table  1,  Figs.  1,2)  which  are  averages  of  all  our  data,  Xe 
was  used  as  the  charge  transfer  gas  in  order  to  observe  the  effects,  if  any, 
of  excited  states  of  oxygen. 


The  attenuation  of  a beam  by  a scattering  gas  is  given  in  the  usual 

I = IQ  exp  (-nQl)  (1) 


notation^  by 


When  the  gas  contains  several  species  with  mole  fraction  xi  we  have  instead 

I = IQ  exp  C-nl  ^ x.?.]  (2) 

which  is  experimentally  indistinguishable  from  the  case  of  a pure  gas. 

(81 

However,  Turner  et  alv  ' have  noted  that  when  the  beam  is  mixed 


1 = lo  $xiexP["nV]- 


(3) 


Here  a plot  of  1 n ( 1/ I ) versus  n is  nonlinear,  and  in  principle  one  may 
obtain  the  individual  by  curve  fitting.  In  a study  of  charge  transfer^ 
the  authors  were  able  to  obtain  cross  sections  for  a mixed  beam  of 
energetic  0+(  S°)  and  0+(20°).  However,  such  an  experiment  is  not  sensitive 
to  multiple  collisions  so  that  90%  or  even  95%  attenuation  could  be  achieved. 
In  the  present  work,  due  to  the  effect  of  multiple  collisions,  we  are 
restricted  to  attenuation  <70%  in  which  range  nonlinearity  is  manifested 
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only  for  a very  large  ratio  of  cross  sections.  It  is  thus  not  surprising 
that  nonlinearity  was  not  observed  in  the  present  experiments.  However, 
provided  the  cross  sections  differ  for  different  states  of  the  beam,  the 
measured  cross  section,  being  a weighted  average  over  states,  will  change 
as  a function  of  relative  concentration.  Careful  measurements  on  the 
scattering  with  Ne  revealed  no  significant  changes  in  cross  section 
irrespective  of  the  method  of  production  of  the  0 atoms,  contrary  to  the 
findings  of  BKL. 

Discussion 

BKL  assumed  that  electron  transfer  from  0"  results  exclusively  in 
0(?o),  the  energy  defect  being  AE  = -1.03eV^  for  the  case 

cf(2P)  + o2(x3e)  - o(3P)  + o2'(x2n).  (4) 

But  we  cannot  dismiss  the  possibility  of  forming  0(*D)  although  since 
AE  = -2.99eV  it  is  less  probable.  In  the  case  of  0+,  charge  transfer  with 
Xe  was  assumed  to  proceed  via 

0+(4S)  + Xe(*S)  - 0(3P)  + Xe+(2P),  AE  = -1.49eV  (5) 

and 

0 (4S)  + Xe( XS)  - 0(!D)  + Xe+(2P),  AE=0.47eV.  (6) 

Charge  transfer  with  Kr  represents  an  intermediate  case  since  the  energy 
defects  for  (5)  and  (6)  are  respectively  0.38  and  2.34eV.  We  note  that 
reaction  (6)  violates  Wigner's  spin  rule.  If  one  invokes  interatomic 
LS  coupling  then  the  observations  of  BKL  are  not  necessarily  inconsistent 
with  ar,  increased  concentration  of  0(3D)  as  one  changes  from  Kr  to  Xe, 
due  to  the  increasing  mass.  However,  this  mechanism  does  not  explain  the 
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independence  of  cross  sections  as  a function  of  charge  transfer  reaction 
in  this  work.  From  the  stated  range  of  the  potential,  determined  by  the 
distance  of  closest  approach  r0,ar,d  the  relation  Q - nro2,  we  deduce  that 
the  O-Ne  cross  sections  of  BKl  changed  by  as  much  as  11%  at  the  lowest 
energy  which  is  we^l  beyond  experimental  error,  but  only  1.4%  at  4000  eV,  whereas 
in  the  present  work  the  maximum  difference  was  only  3.6%  throughout  the  range. 

The  0-Ne  cross  sections  with  Kr  and  Xe  charge  transfer  were  here  an  average 
of  1.6%  and  1.3%  respectively  lower  than  for  0"  with  02>  but  this  is 
within  experimental  error  estimated  as  + 1%.  Indeed  scattering  with 
a He  target  gas  produced  an  opposite  but  equally  insignificant  trend.  As 
a further  check  a number  of  measurements  were  made  with  N2  and  also  H^O  as 
charge  exchange  gas  yielding  the  same  cross  sections.  The  energy  defects 
for  production  of  0(3P)  and  O^D)  are  respectively  1.97  and  3.93 eV  for 
N2>  and  -1.00  and  + 0.97eV  in  the  case  of  HgO.  Again  the  reactions 

0+  (4S)  + N2  (X1!)  - 0(XD)  + N2+  (X2£)  (7) 

and 

0+  (4S)  + H?0  (X^Ai)  - 0( lD)  + H20+  (8) 

are  spin  forbidden.  Charge  transfer  reactions  are  in  general  highly 
velocity  dependent,  but  five  differert  reactions  yielded  the  same  0-Ne 
cross  sections  over  the  entire  energy  range.  Theoretical  calculations^^ 
give  a much  higher  potential  for  the  triplet  state.  We  thus  conclude 
that  our  beam  contained  a negligible  fraction  of  the  metastable  state.. 

The  results  provide  strong  evidence  that  spin  is  conserved  in  these 
experiments,  and  also  then  in  those  of  BKL  since  their  lab.  energy 
extending  from  290-4000eV*  largely  overlaps  this  work.  We  must  therefore 
seek  an  alternative  explanation  of  their  results.  One  possibility  is  that 


* It  appears  that  the  stated  lowest  energy  of .r2|°v\\^/^&na%0tfower 
system.  From  the  relation  InQ/lnE-  -2/s  where  V Kr 
limit  of  600eV  for  the  0-Ne  systems  of  BKL.  The  range  600-4000 eV  was 
normally  used  in  previous  work  (Ref. 11).  IV-5 


their  electron  multiplier^)  is  more  sensitive  to  metastables  than  the 
Channel tron  used  here,  which  would  explain  the  increasing  effect  with 
lower  energy.  However,  the  excitation  energy  of  0(  lD)  is  only  1.96eV 
which  seems  insufficient  for  the  production  of  Auger  electrons.  The  major 
difference  in  the  cross  sections  of  BKL  is  between  0+  and  0"  rather  than 
the  different  charge  transfer  reactions  of  0+.  Such  differences  may  arise 
if  the  ion  extraction  field  penetrates  the  source  plasma  to  the  cathode 
region  resulting  in  a higher  energy  for  0 and  a lower  energy  for  0+. 

This  is  more  significant  at  lower  energies,  but  can  be  dismissed  since  the 
predicted  deviations  are  the  opposite  of  those  observed.  A more  probable 
cause  also  lies  in  the  different  design  and  operation  of  the  ion  source. 

We  have  already  remarked  that  Turner  et  al  found  up  to  30%  0+(2D)  as  the 
electron  energy  was  increased  from  threshold  up  to  lOOeV.  Due  to  the 
different  extraction  efficiencies  for  the  and  2D  states,  resulting  from 
formation  on  different  repulsive  states  of  02+,  the  relative  concentrations 
may  be  a sensitive  function  of  the  source  construction. ^ Charge  transfer 
of  0+  (2d)  with  Kr,  Xe,  N2  or  H20  may  yield  either  0(3P)  or  O^D)  with 
conservation  of  spin.  The  energy  defects  are  respectively  -2.92,  -4.79, 
-1.34  and  4.30eV  to  form  0(3P),  and  -0.96,  -2.83,  +0.62  and  -2.34eV  to 
give  0(1D).  Thus  excited  0+  ions  may  preferentially  form  0( 1 D) , although  in 
the  case  of  N2  we  might  also  expect  to  form  the  excite<j  N'2  (A2nu)  or 
vibrational ly  excited  n£  (XZE)  which  favor  0(3P).  Similarly,  H20+ 

( A2 A j ) may  participated  In  this  work  we  were  unable  to  detect  any  changes 


T We  may  discount  excited  states  of  Kr+  or  Xe+  since  they  lie  c.a.  9eV 
above  the  ground  state. 
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in  cross  section  when  source  operating  conditions  were  altered.  The 
variable  parameters  were  source  pressure,  electron  energy,  extraction 
potential,  and  the  use  of  CO  instead  of  Og  in  the  source.  This  may  simply 
imply  that  excited  0 ions  were  present  in  this  work,  but  that  despite  the 
unfavorable  energy  defect  the  dominant  channel  is  formation  of  0(^P), 
making  them  indistinguishable  from  the  ground  state  ions.  The  explanation 
of  BKLs ' observations  is  thus  not  definitely  established,  but  may  be  due  to 
the  presence  of  0+  rD)  in  their  ion  beam. 

Although  we  have  assumed  that  the  relative  differences  in  the  O-Ne 
cross  sections  of  BKL  are  beyond  experimental  error,  such  differences  in 
absolute  determinations  by  independent  workers  are  not  uncommon.  This  is 
evident  from  the  difference  between  BKL  and  this  work  in  the  O-Ne  potential 
obtained  from  0 (Fig. 2)  where  there  can  be  no  difference  in  the  relative 
proportions  of  singlet  and  triplet  atoms  at  a given  lab.  energy.* 

Leonas,  in  a recent  review,'  has  corrected  the  earlier  measurements  of 
all  the  oxygen-noble  gas  systems^’^’  for  errors  in  the  size  of  the 
aperture,  which  changes  the  magnitude  of  the  potential  but  leaves  the  range 
unchanged.  Agreement  with  this  work,  however,  is  not  improved  (Figs.  1,2). 


Assuming  that  the  relative  concentrations  of  0(^P)  and  0(!D)  change 
as  a function  of  lab.  energy,  a difference  between  the  potentials  of  BKL 
and  this  work  may  result  since,  due  to  differences  in  angular  resolution, 
the  same  lab.  energy  gives  information  about  different  regions  of  the 
potential  curve  in  the  two  experiments.  However,  the  reported  differences 
in  potential  seem  too  large  to  be  accounted  for  in  this  way. 
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Leonas  notes^’^  that  calculations  of  the  molecular  singlet  state  by 

Allen  et  al.^^  lie  substantially  above  the  measurements  of  BKL.  Since 

the  triplet  lies  at  still  higher  energies,  this  is  not  consistent  with  the 

view  that  the  measurements  represent  some  mixture  of  states.  He  thus 

concludes  that  the  calculations  are  in  error.  The  present  work  suggests 

that  the  theory  is  more  likely  to  be  correct.  The  authors^)  also 

estimated  that  the  molecular  triplets  lie  even  higher  than  the  doublets  of 

the  corresponding  fluorides.  The  comparison  with  this  work  for  Ne  and  He 

(Figs.  1,2)  shows  that  they  are  nearly  equal  which  is  reasonable  agreement.* 

A further  estimate  of  the  triplet  states  is  provided  by  Thomas- Fermi -Dirac 

(TFD)  model  calculations. The  heteronuclear  diatomirs  are  obtained 

from  the  corresponding  homonuclear  potentials  by  the  empirical  combination 

1/2 

rule  V12  - (vn V2?)  , We  see  (Figs.  1,2)  that  the  results  of  Gaydaenko 

and  Nikulin  are  in  excellent  agreement  with  this  work,  whereas  the 
calculations  of  Abrahamson  are  known  to  be  too  high .06,17) 


* 

The  unpublished  calculations  of  the  triplet  states  must  be  regarded 
with  caution  since  they  were  at  the  limit  of  technical  capabilities 
then  available,  but  in  this  range  they  certainly  lie  above  the 
singlet  state  calculations,  which  are  considered  to  be  reliable. 
L.C.  Allen  (private  communication). 
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The  O-Kr  and  O-Xe  potentials  of  Belyaev  and  Leonas^)  are  also 
corrected  for  errors  in  pressure  calibration,  which  does  alter  the  range.* 
Comparing  these  revised  potentials  (Fig. 3)  we  see  that  they  all  lie 
substantially  below  those  obtained  in  this  work,  but  the  difference  is 
probably  too  large  to  be  explained  on  the  basis  of  excited  0 atoms. 
Furthermore,  instead  of  the  expected  increase  in  potential  with  size  of 
the  target  atom,  they  show  reverse  ordering  for  Ar,  Kr  and  Xe,  and  predict 
the  same  potential  for  O-Ne  and  O-Ar  at  a sep*  ation  of  1.88.  Possibly, 
at  least  in  the  case  of  Xe,  there  is  typographical  error;  increasing 
: the  O-Xe  potential  by  a factor  of  10  brings  it  into  moderate  agreement 

with  this  work.  The  TFD  calculations  of  Gaydaenko  and  Nikulin  are  in  fair 
H agreement  with  this  work. 


I-  I 
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does  not  qVote  the  ran9e  of  validity.  We  must  therefore  assume 
greatlyealteredSeCt1°n  ^ thUS  th*  range’  which  varies  as  is  not 
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Summary 


We  find  that  the  oxygen  beam  used  in  this  work  consists  only  of 
ground  state  atoms,  irrespective  of  the  charge  transfer  reaction  used  to 
generate  that  beam.  This  is  still  consistent  with  the  observations  by 
BKL  of  a mixture  of  ground  and  excited  states  from  the  same  reactions. 

However,  the  formation  of  excited  0(1D)  atoms  from  ground  state  ions  is 
spin-forbidden  and  does  not  occur  to  any  significant  extent.  We  postulate 
that  the  ion  beam  of  BKL  contained  excited  0+(^D)  ions.  In  future  experiments 
a pure  0(  P)  beam  might  best  be  obtained  by  photodetachment  from  C~.  It 
has  been  suggested^)  that  this  method  might  also  be  used  to  obtain 
fast  neutral  beams  with  a high  concentration  of  *D  and  *S  states  for  which 
the  threshold  energies  are  3.433eV  (361 l8)  and  5.656eV(21928)  compared 
with  1.465eV( 84638)  for  the  ground  state.  By  calibration  against  a reliable 
theoretical  potential  for  say  0(lD)-He^^  one  could  then  obtain  the 
interactions  with  0(1D)  for  any  other  system. 

The  resulting  oxygen-noble  gas  potentials  are  much  larger  than  those 
of  BKL  and  still  higher  than  the  revised  estimates  of  Leonas.  The  absolute 
difference  seems  too  great  to  be  attributed  solely  to  excited  states, 
but  theoretical  estimates  are  in  better  agreement  with  the  present  work. 


II  Atmospheric  Molecules 


Analysis  and  results 

The  interactions  of  molecules  are  characterized  by  an  anisotropic 
potential,  which  may,  however,  in  the  simplest  approximation  be  treated  as 
a spherically  symmetric  potential.  The  analysis  which  then  proceeds 
exactly  as  for  atoms  will  be  referred  to  here  as  the  potential  averaging 
(PA)  procedure.  Until  the  recent  work  of  Leonas^  nearly  all  scattering 
experiments  were  so  treated. Theoretical  calculations  on  the 
He-I^  ^ and  systems  have  used  a one-center  expansion  in  Legendre 

polynomials  to  describe  the  anisotropic  potential: 

V(R,d,tii)  =£  V^(R,d)  P^ ( cos  ij<)  (9) 

where  R is  the  distance  between  centers  of  mass  of  the  atom  and  molecule, 
d is  the  molecular  bond  length  and  ip  is  the  angle  between  R and  d.  Due 
to  the  difficulty  of  uniquely  determining  more  than  two  parameters, 
experimental  work  has  utilized  the  additive  potential  model  in  which  the 
intermolecular  potential  is  simply  the  sum  of  the  individual  interatomic 
interactions.  Thus  for  the  system  A + BC  we  have 

V(R.^)=  vAB(rAB)  + vftc(rAC)  (10) 

for  a particular  orientation,  and  the  spherically  symmetric  potential  \T( rJ 
is  the  average  of  all  orientations  of  V(R^).  Amdur,  Mason  and  Jordar/21  *22) 
have  successfully  used  this  model  to  extract  effective  atom-atom  interactions 
of  the  form  V = Kr  s (where  K and  s are  constants)  from  experimental  PA 
interactions  of  atoms  with  homonuclear  diatomics.  With  the  aid  of  combination 
rules  they  then  construct  molecule-molecule  potentials  in  good  agreement 
with  other  data.  The  extension  to  heteronuclear  diatomics  and  more 
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complicated  systems  requires,  in  general,  the  determination  of  more  than 

two  potential  parameters  from  the  experimental  data.  Leonas  overcomes 

this  difficulty  by  assuming  that  all  atoms  within  a molecule  are  identical. 
(23) 

Amdur  et  al ' 7 have  introduced  the  peripheral  force  approximation  to 

treat  CH^,  in  which  the  interaction  of  the  central  carbon  atom  is  ignored. 

Sinai.oglu  has  noted^  that  in  the  time  scale  of  a collision,  the 
target  molecule  undergoes  negligible  rotation  under  these  experimental 
conditions.  Given  a trial  potential,  one  should  therefore  calculate  the 
cross  section  for  each  orientation  and  then  average  to  compare  with 
experiment.  This  procedure  is  termed  cross  section  averaging^8) (CSA)  to  dis- 
tinguish it  from  the  straightforward  PA  method.  The  latter  would  be 
valid  in  the  limit  of  rapid  rotation  of  the  target.  Expressing  the 
potential  as  an  expansion  in  Gegenbauer  (ultra-spherical)  polynomials*, 
Sinano^lu  estimates^  that  the  Ar-N2  potential  of  Amdur^2^  is  increased 
by  about  5%  in  this  way,  which  is  within  experimental  error.  The  present 
work  extends  to  shorter  range  where  the  difference  may  become  larger 
as  the  molecular  structure  becomes  increasingly  important.  Furthermore, 

Leonas  states  that  the  two  methods  differ  by  a factor  of  two  in  some 
(5) 

cases  which  is  certainly  well  beyond  experimental  error. 

We  shall  treat  only  the  interactions  of  0 atoms  with  H^,  and  0^ 
by  the  CSA  procedure.  In  order  that  other  methods  in  addition  to  PA 
may  be  applied  to  these  and  the  remaining  data,  the  cross  sections  ^xp(E), 
including  an  end  correction,  are  given  in  Table  II.  The  scattering  geometry 
is  given  above.  PA  parameters  A and  a are  given  in  Table  III  and  the  results 
compared  with  other  data  in  Figs.  4-9.  The  only  other  direct  experimental 


* 

This  is  equivalent  to  the  assumption  of  additive  inverse  power 
interatomic  potentials.  ^2^ 
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Their  PA  curves  lie 


data  are  those  of  Leonas  and  coworkers. ( 5 ,25-29 ) 
both  above  and  below  this  work,  the  largest  disagreement  being  50- 1 00% 
in  the  cases  of  Ng  and  (Figs. 5, 6).  Jordan,  et  al^22^  indirectly 
predict  potentials  for  0-^,  Og  and  NO  by  applying  an  additive  potential 
model  together  with  the  geometric  mean  combination  rule  to  their  PA 
results  for  Ar-Ar,  N£  and  Og.  These  are  slightly  higher  than  this  work 
but  in  reasonable  agreement,  whereas  the  semi -empirical  calculations  of 
Yun  and  Mason, ^ and  of  Vanderslice  et  al^31^  are  probably  too  high  (Figs. 5-7). 

Considering  now  the  CSA  method,  we  assume  that  the  interatomic 
potential  for  say  0-H  in  the  O-H^  system  is  given  by  V(r)  =A  exp(-ar). 

In  addition  to  averaging  over  the  scattering  path  at  each  energy,  as  in 
the  central  field  or  PA  case,  we  must  now  integrate  over  all  orientations 
of  the  target  to  obtain  the  calculated  cross  section  (Jca^c(E).  A two 
parameter  search  for  A and  a would  be  a costly  computation.  However,  a 
simple  iterative  procedure  converges  rapidly  to  the  one-parameter  search 
for  A independent  of  a as  in  the  case  of  the  noble  gases  (see  Appendix  B). 

The  spherically  averaged  atom-molecule  potential  is  then  given  by^32^ 

VUTT  = 2Ae'0tR(a2Rd)'1  [2(aR+l)  sinh(ad/2)  -ad  cosh  (ad/2)].  W) 

In  practice  this  function  may  be  approximated  by  the  form 

V(R)  = A exp  (-a  R) 


★ 

From  experimental  considerations,  the  Ar-N„  potential  in  Ref.  22,  and 
hence  those  derived  from  it,  are  known  to  be  somewhat  too  high. 

J.E.  Jordan  (private  comnunication) . 
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with  negligible  error.  We  take  tse  range  of  validity  to 
be  the  same  as  that  given  by  potential  averaging.  The  range  of  the 
"interatomic"  potentials  V ( r)  =A  exp(-ar)  is  not  obvious.  However,  we 
follow  Vanderslice  et  al  and  substitute  the  limits  of  R in 

n _ Pi  , / d\2i  Mr 


R = r [1  + TT 


to  obtain  the  limits  of  r.  The  results  are  given  in  Tables  IV  & V. 


Discussion 

The  CSA  potentials  for  H2,  N2  and  G2  (Figs.  4-6)  are  slightly  higher 
than  the  PA  results  and  in  excellent  agreement  with  the  predictions  of 
Sinano^lu.  We  find  that  the  differences  increase  with  increasing  bond 
length  and  decreasing  R as  expected.  However,  even  for  02  the  maximum 
difference  is  only  18£.  The  foregoing  comments  apply  equally  well  to 
the  "interatomic"  potentials  (Figs.  10-12)  which  comprise  the  final 
intermolecular  potentials.  Here  the  PA  pofentials  were  obtained  by  a 
method  analagous  to  that  employed  by  Amdur  et  al  for  an  inverse  power 


representation. 


The  parameters  are  listed  in  Table  VI.  Clearly 


these  potentials  cannot  be  expected  to  apply  to  free  atoms,  but  we 
nevertheless  present  them  for  comparison  with  similar  data.  As  required 
for  internal  consistency,  the  results  of  Jordan  are  again  slightly  higher, 
and  the  semi-empin'cal  calculations  still  larger,  with  the  exception  of 

the  0-N  reSUlt  of  Meador  1 annac'  nntontialc  ann  murh  lnyor  than  the  » 


Leonas'  potentials  are  much  lower  than  the  rest. 


We  are  unable  to  explain  the  large  differences  found  by  Leonas 
between  the  PA  and  CSA  procedures.  In  particular,  his  finding  of  lower 
CSA  than  PA  potentials  (for  CO  they  intersect  and  for  C02  the  earlier 
CSA  results  of  Kalinin,  et  ai.(28,29)  are  higher)  is  contrary  to  the 


IV- 1 4 


predictions  of  Sinanotflu  and  results  of  this  work.  Since  the  CSA  results 

appeared  in  a later  review^  it  is  possible  that  changes  in  calibration 

altered  the  previous  experimental  cross  sections.  Comparison  is  also 

hindered  by  a failure  to  state  the  range  of  validity.  It  is  implied  that 

the  potentials  are  valid  in  the  range  0.1  - lOev/5,29,33^  but  assuming 

that  the  CSA  results  are  a reanalysis  of  the  earlier  PA  data,★ **  then  this 
must  overstate  the  range.  We  conclude  from  a comparison  of  this  and  other 

i i 

work,  that  at  least  in  the  case  of  O-O2  and  0-^,  the  CSA  potentials  of 
Leonas  are  incorrect. 

The  usefulness  of  the  additive  force  model  is  demonstrated  by  the 
0-N0  potential.  Combining  our  CSA  results  for  0-0  and  0-N  obtained  from 
O-O2  and  O-N2,  we  predict  a potential  for  0-N0  which  is  in  very  close 
agreement  with  our  direct  PA  result  (Fig. 7).  This  method  is  more  accurate 
than  the  alternative  assumption  that  the  molecules  have  spherical  symmetry. 

This  is  demonstrated  in  Fig. 8 where  we  show  an  indirect  prediction  of  the 
0-C0  potential  (V(R)  = 539  exp(-2.932R)eV,l .31sR<2.448)  obtained  by 
applying  the  geometrical  mean  combination  rule  directly  to  the  potential 
averaged  Ar-CO,  Ar-Ar  and  O-Ar  interactions.*** 


★ 

Changes  in  calibration  of  geometry  and  pressure  were  applied  to  the 
0-noble  gas  data. 

** 

Such  an  assumption  is  made  on  the  grounds  that  these  and  all  the  other 
data  (with  the  apparent  exception  of  0-Hp  which  we  have  not  found  in  the 
literature)  have  been  previously  published  in  PA  form. 

*** 

The  Ar-Ar  potential  and  the  potential  averaged  Ar-CO  interaction  obtained 
in  this  laboratory  are  respectively  V(r)  = 5061  exp(-3.598r)eV, 

1.70  sr  s 2.65A  and  V(R)  = 3450  exp(-3.313R)eV,l .79  ^ R c 2.76A. 
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We  have  not  thought  it  worthwhile  to  make  the  more  drastic  assumptions 
required  for  a direct  CSA  treatment  of  the  heteronuclear  diatomics  and 
triatomics.  Clearly  as  the  molecule  becomes  larger  and  more  complex 
the  differences  between  PA  and  CSA  potentials  become  more  severe. 

However,  with  the  caution  that  they  are  probably  at  least  10%  too  low, 
the  PA  potentials  obtained  here  for  0 atoms  interacting  with  NO,  CO, 

C02  and  H20  will  be  adequate  for  many  purposes.  We  note  that  the  difference 
between  the  PA  potentials  is  small  for  the  isoelectronic  0-N2  and 
0-C0  interactions. 

One  must  also  remember  that  chemical  interaction  may  be  invol/ed 

in  all  the  oxygen  atom-molecule  collisions  so  that  the  potentials  are  also 

an  average  over  states.  The  higher  angular  resolution  employed  by  Leonas 

(5  25) 

has  permitted  the  observation  of  such  curve  crossings.  ’ 


APPENDIX  A 


Analysis  of  atom-atom  scattering 

For  systems  involving  a projectile  of  mass  m-|  and  target  of  mass  m2, 
where  m-j  < m2,  the  traditional  analysis  in  terms  of  the  first  order 
momentum  approximation^ is  entirely  adequate.  However,  when 
m-|  > m2  (e.g.  O-He  or  the  potential  averaged  0-H2)  it  becomes  necessary 
to  include  higher  order  terms  in  the  momentum  expansion.  These  terms 
have  been  evaluated  by  Lehmann  and  Leibfried^^^  and  by  Smith  et  al^) 
for  several  forms  of  the  potential  V(r).  It,  an  exact  analysis  the 
classical  center  of  mass  (c.m.)  scattering  angle  x is  related  to  the 
impact  parameter  b by 

* ‘ *-2b/r  t1  -$r  7 (a” 

o 

where  rQ  is  the  largest  zero  of 

VM  * 4 - 1 (A2) 

tT 

and  E-j-  is  the  total  energy  of  the  system.  Here  we  substitute  for  ET  the 
initial  kinetic  energy  Er  in  the  absence  of  interaction  given  by 

Er  = m2  E/ (m-j  + m2)  (A3) 

where  E is  the  projectile  lab  energy.  In  these  experiments  the  target 
thermal  energy  is  negligible. 

The  momentum  approximation  sets  b = rQ  (valid  when  V « Ep)  and  one 
evaluates  Eq.  A1  as  an  expansion  in  powers  of  V/Ep: 


X = s X-i 
i 


For  the  exponential  form  V(r)  = A exp(-ar)  the  first  few  terms  are 


X1  = ab  KQ(ab) 


= - — K,  (2ab) 
Er 


X3  = | (ab)3  K0(3ab)  - Sab^Oab)  + ^ (ab)2  KQ(3ab)  (A7) 

Er 


where  is  the  ith  order  modified  Bessel  function.  In  the  usual 


i 

analysis^3, 34^  6 , the  minimum  lab  angle  required  for  the  projectile  to 


miss  the  detector  after  scattering  from  a given  point  in  the  cell,  is 


related  to  x,^  via  the  small  angle  approximation 


m-|  + mg 


The  corresponding  impact  parameter  bm  is  then  obtained  from  Eq.A5  by 


Newton-Raphson  iteration  to  give  the  single  ray  cross  section  (ignoring 
back  scattering) 


qca1c(E)  = irCbjE,^)]2 


for  that  point.  The  final  calculated  cross  section  Q (E)  is  then  the 


volume  average  of  Qca1c(E)  over  the  scattering  path.  This  is  evaluated 


ty  a 3-fold  Gauss-Legendre  quadrature,  each  of  8 points  for  a total  of 


512.  The  parameters  A and  a must  be  adjusted  to  give  the  closest  fit  tc 


the  experimental  energy  dependence  of  the  cross  section.  Clearly  as  one 


I 
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includes  higher  order  terms  in  the  expansion  of  x>  the  computation  becomes 
very  long.  Furthermore  the  asymptotic  expression 


>'  ■©'  *'  U • f "$) 

is  no  longer  useful  to  predict  the  value  of  a when  mg  » m-j . 
One  must  then  also  use  the  exact  c.m.-lab  transformation 


(A10) 


tan  0 = 


m«  sin  x 


(All) 


mg  cos  x + mi 

For  the  present  geometry  the  maximum  lab  angle  detected  is  0.0285  rad. 

Xe  - He  at  100  eV,  to  consider  an  extreme  case,  then  giver  the  corresponding 
forward  c.m.  angle  x = 1*24  rads,  which  is  no  longer  small.  Use  of  Eq.A8 
gives  22%  error  in  x*  This  could  be  reduced  to  1.6%  by  using  the  second 
order  expression  for  forward  scattering 


X = arcs in 


-cot  0 + cot  0 1 + 2 


2 1 ^ 
tan  0 


(A12) 


However,  in  this  work  Eq.  All  was  solved  iteratively  for  x starting 
from  (A12). 

Analysis  for  an  inverse  power  potential  is  much  more  direct  giving 
to  first  order 

/Kf  \2/$ 

QCa1c(E)  = < 0m'2/s  > (A13) 

= Ji 


where  C1  - ir  r(s/2  + 1/2)/  r(s/2)  and  to  2nd  order 


qcalC(E) 


fa}\2/s 


rn 


} 


<e. 


m 


-2/s.  2^ml  + ^ 
- 2 


sni2C  i 


c9  <0  1_2/s> 

i m 


(A14) 
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where  C, 


(3) 


ir  r(s  + %)/2r(s  - l).v  ; However,  (A14)  assumes  the  validity 


of  (A8)  which  is  inconsistent  with  the  need  to  include  X2  in  the 


momentum 


expansion.  It  would  be  preferable  to  use  x(0m)  from  (A1 2)  in  the  expression 


qca,c(E 


>•[ 


(m1  + m2)KC1-j2/s 


m2E 


] { 


<[x(em)]-2/s>-^-2<[x(8m)]1-2/s>; 

SL1 


(A15) 


The  convergence  of  Eq.  A4  has  been  considered  in  detail  by  Leibfried 
(37) 


and  Plesser.w''  We  see  from  Fig.  13  that  with  this  apparatus  the  series 

diverged  at  all  energies  for  Xe  - He  (m2/m.|  - 32).  Even  the  Ar  - He 

system  (m^/m^  = 10)  diverges  here  at  the  lowest  energies:  this  case  has 

formerly  been  analyzed  in  terms  of  Eq.  A13.^)  Thus  caution  must  be 

exercised  in  using  the  momentum  approximation  particularly  since  the 

leading  term  gives  results  which  may  appear  reasonable  in  the  absence  of 

3/2 

data  for  comparison  (see  Fig.  14).  Furthermore,  for  ab  > 3 the  series 
only  converges  in  alternating  groups  so  that  examination  of  the  first  few 
terms  is  not  a sufficient  test  of  convergence. 

When  the  series  (A4)  diverges  it  is  necessary  to  evaluate  Eq.  A1 
exactly.  In  the  case  of  potential  averaged  O-H^  (m2/m-|  - 8)  it  is 
desirable  since  convergence  is  slow.  This  is  conveniently  achieved  by 


setting  y = r/r  which  transforms  (Al)  to  the  interval  [0,1]  in  y,  and 


(39) 


using  an  8-point  Gauss-Tchebyshev  quadrature'  or  alternatively  Gauss- 


Mehler  quadrature to  remove  the  square  root  singularity.  In  some 

(3) 


cases  we  have  used  the  more  rapidly  convergent  form  of  Eq.  Al 


X = 


c 

>/ 


r dV 
(Ej-V ) dr 


H-S] 


dr 

T 

r 


(A16) 
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We  require  the  value  of  b given  x.  but  the  Newton-Raphson  method  is  no 
longer  useful.  Instead  we  set  up  a table  of  log  x versus  b for  a given 

ratio  A/Er  and  obtain  b by  inverse  interpolation,  20  or  perhaps  40  points 

being  adequate  in  most  cases. 

We  may  reduce  the  two  parameter  search  to  one  by  introducing  reduced 
variables.  For  the  exponential  form  we  have 

p = ar 

3 = ab  (A17) 

o = a2Q 


and  V(p)  = Ae~p.*  The  dependence  on  a of  Eqs.  A1  and  A15  is  thus  eliminated 


In  practice  we  therefore  choose  a value  of  A and  for  each  energy  calculate 
the  scattering  path  volume  average  5ca1c(Er)  of  the  single  ray  reduced 
cross  sections  oca1c(Er)  = The  reduced  impact  parameter  is 

evaluated  by  inverse  interpolation  of  log  x versus  6 from  (A18).  a is 
then  given  by  the  least-squares  condition 


* The  exponential  potential  may  be  unrealistic  for  small  molecules  at  very 
short  range  (e.g.  H-He)  so  we  have  followed  Jordan  and  Amdur  in  matching 
the  exponential  to  a screened  Coulomb  potential  at  small  p (<  0.6). 
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- calc  - calc 

o„  o. 


tr.  vt 

i ;*ca,c 


(A19) 


where  the  suffix  l denotes  summation  over  all  energies.  A is  adjusted  to 
minimize  the  sum  of  squares  of  deviations  x (not  to  be  confused  with  the 
c.m.  scattering  angle)  between  experimental  and  calculated  cross  sections. 
Once  three  trials  of  A have  been  executed,  a rapid  convergence  is  obtained 
by  assuming  (A.,  xf)  satisfy  a quadratic  in  A,  the  minimum  giving  a new 

prediction  of  A. 

This  technique  is  equally  applicable  to  the  momentum  approximation, 
but  due  to  the  increased  efficiency  the  latter  offers  no  advantages  over 
an  exact  evaluation  of  Eq.  A1  even  when  m1  4 m2  and  the  convergence  of 
Eq.  A4  is  good.  A program  written  in  Fortran  IV  performs  the  exact 
calculation  in  0.080  pq  secs,  of  IBM  370/65  central  processor  time, 
where  p and  q are  respectively  the  number  of  energies  and  interpolation 
points.  Thus  the  complete  determination  of  a potential 
requires  typically  30  seconds  of  computation. 

Applying  this  method  to  He-Xe  gives  indistinguishable  results  from 
the  1st  order  momentum  approximation  (Fig.  14).  However,  for  Xe-He  there 
is  a marked  improvement  in  the  exact  result.  The  agreement  with  the  He-Xe 

result  of  Amdur  end  Engler'43'  and  the  TFD  calculation  of  Abrahmson04' 

1,  also  good.  In  the  cases  of  0-H2  and  O-He,  the  changes  are  significant 

but  less  dramatic. 
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APPENDIX  B 

Cross  Section  Averaging 

We  assume  that  the  i ntermol ecul ar  potential  may  be 
represented  by  the  sum  of  the  interatomic  teams  given  in 
exponential  form  V = Aexp(-p^)  (see  Fig.  15),  thus 

V( p.uO  = A ( e P'  + e ^2)  (Bl) 


where 


Pl  = (p2  + 62  - 2p6cosij») 


p p 

= (p  + 6 + 2p6cos\J> ) 


For  a given  energy  and  orientation  of  \p  we  may  substitute 
V(Pti|>)  into  Eq.  A17  and  proceed  as  before  to  calculate 
the  volume  average  oca^c(^).  Averaging  over  orientations 
we  then  obtain 

5calc  = /j5ca'c(*)d(cos*).  (B3 ) 

Here  we  use  Gauss-legendre  quadrature  giving  an  8-fold 

increase  in  computation  time  over  the  atom -atom  case. 

Since  Eqs.  B2  require  the  reduced  rather  than  the  actual 

known  bond  length  we  cannot  eliminate  the  interatomic 

potential  parameter  a as  before.  However,  given  a trial 

value  ai  we  may  obtain  a new  estimate  using  the  condition 

(A19),  namely  1/2 

EoCSlC(at )oCa1C (a. ) 

-,L± 5 . (B4) 

EQcxP5calc(a.) 


, 
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One  or  occasionally  two  iterations  are  usually  sufficient 
to  give  4 sig.  figs,  in  a,  so  that  the  problem  is  effectively 
reduced  to  a one-parameter  search  for  A.  Preliminary  guesses 
at  A and  a may,  if  desired,  be  obtained  from  the  PA  method 
(see  main  text).  Typically  4 mins,  are  required  for  an  IBM 
370/65  computer  to  calculate  the  optimum  A and  a in  these 
experiments.  The  atom-molecule  potential  is  given  by  a 
spherical  average  of  the  summed  interatomic  contributions: 

TO  = 2Ae"aR(a2Rd)"1[2(aR+l)sinh(ctd/2)-adcosh(ad/2)].  (B5) 

For  the  systems  studied  here,  V (R)  may  be  closely  represented 
by  an  exponential  form  V ( R ) = Aexp(-aR).  The  range  of 
validity  is  taken  to  be  the  same  as  that  given  in  the  straight- 
forward potential  averaging  analysis. 
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Table  I.  Potential  parameters  for  interaction  of  0(  P)  with  noble  gases, 
V (r ) =A  exp  (-ar) 


Target 

A (eV ) 

<«-') 

Range  (A) 

He 

378 

3.744 

0.85  - 1.81 

Ne 

1222 

3.917 

1.16  - 2.05 

Ar 

791 

3.217 

1.25  - 2.36 

Kr 

1063 

3.249 

1.35  - 2.43 

Xe 

1432 

3.287 

1.43  - 2.50 

Table  II.  Dependence  of  incomplete  integral  cross  sections  (including  end 
correction)  on  lab.  energy  for  collisions  of  0(JP)  with 
atmospheric  molecules. 


Energy(eV) 

Cross 

section  (A2) 

H2 

N2 

°2 

NO 

CO 

co2 

h2o 

2500 

2.25 

4.76 

5.54 

5.34 

4.98 

7.49 

5.01 

2000 

2.56 

5.28 

6.23 

5.89 

5.48 

8.31 

5.50 

1500 

2.99 

5.98 

7.06 

6.62 

6.12 

9.26 

6.14 

1000 

3.65 

7.20 

8.38 

7.94 

7.46 

11.05 

7.40 

700 

4.56 

8.51 

9.68 

9.37 

8.76 

12.85 

8.74 

500 

4.97 

9.92 

10.87 

10.43 

9.96 

14.47 

9.67 

400 

5.49 

10.94 

11.54 

11.31 

10.89 

15.43 

10.41 

300 

6.21 

12.35 

12.70 

12.52 

12.13 

16.74 

11.28 

200 

7.22 

14.29 

14.13 

14.41 

13.58 

18.32 

12.71 

150 

8.11 

15.74 

15.04 

15.50 

14.94 

19.43 

13.78 

100 

9.60 

17.62 

16.25 

16.98 

16.77 

20.97 

15.02 

Table  III.  Potential  averaged  (PA)  parameters  for  interaction  of  0(^p)  with 
atmospheric  molecules.  (The  latter  assumed  to  have  spherical 
symmetry).  V(R)  =A  exp  (-aR) 


Target 

A(ev) 

a(M) 

Range  (A) 

H2 

422 

3.912 

0.84  - 

1.76 

"2 

541 

3.012 

1.20  •• 

2.38 

°2 

1870 

3.670 

1.36  - 

2.31 

NO 

1048 

3.350 

1.30  - 

2.35 

CO 

825 

3.286 

1.24  - 

2.32 

co2 

2347 

3.317 

1.59  - 

2.63 

HP0 

1484 

3.724 

1.28  - 

2.21 

Table  IV.  Cross  section  averaged  (CSA)  parameters  for  interatomic  potentials 
V(r)  =A  exp(-ar)  in  collisions  of  0(JP)  with  homonuclear 
diatomics.  Additive  forces  assumed. 


Target 

Interaction 

A(eV) 

Range  (A) 

H2 

0....H 

209 

4.030 

0.90  - 1.78 

n2 

0....N 

316 

3.221 

1.28  - 2.42 

o. 

0 0 

905 

3.856 

1.45  - 2.37 
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Table  V.  CSA  parameters  for  interaction  V( R)  = Aexp(-aR)  of  0(  P)  with 
diatomic  molecules.  Parameters  obtained  by  averaging  over  all 
orientations  the  interatomic  potentials  in  Table  IV. 

Additive  forces  assumed. 


Target 

A (eV) 

a (A-1) 

Range  (A) 

H2 

444 

3.916 

0.84  - 1.76 

N 

698 

3.109 

1.20  - 2.38 

° 

2438 

3.715 

1.36  - 2.31 

NO^ 

972 

3.326 

1.45  - 2.44 

(a)  Note  that  this  potential  is  not  derived  from  0-N0 
experimental  data,  but  rather  from  0-0^  and  0-N,,. 


Table  VI.  Interatomic  PA  parameters  for  V(r)  =A  exp(-ar),  derived  from 
0(3p)  - homonuclear  diatomic  potentials  in  Table  III. 
Additive  forces  assumed. 


Target 

Interaction 

A(eV) 

a(A_1) 

Range  (A) 

H2 

0....H 

202 

4.037 

0.90  - 1.78 

n 

0....N 

249 

3.121 

1.28  - 2.42 

O 

0 0 

703 

3.810 

1.45  - 2.37 
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Comparison  of  experimental'  potentials  for  O-Ar,  Kr,  Xe  with  calculations 
in  the  Thomas-Fermi-Dirac  approx.  0 4,15)  j^e  earlier  results  of 
Belyaev  et  al . (12,13)  are  omitted  to  avoid  confusion;  they  lie  between  this 
work  and  the  results  of  '.eonas. 
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CENTER  OF  MASS  SEPARATION  (A) 


4.  Comparison  of  experimental  potentials  for  0-H?  treated  in  the  potential 
average  (PA)  and  cross  section  average  (CSA)  approx. 
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derived  from  Ar-(L,  Ar-N 

and  Ar-Ar  assuming  additive  interatomic  forces  and  the  geometric^ 
mean  (GM)  combination  rule. 
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interactions  assuming  additive  interatomic  forces, 

.(22)  additionally  assumes  the  GM  rule. 
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CENTER  OF  MASS  SEPARATION  (A) 


Comparison  of  direct  PA  and  CSA  exptl.  results(5’26J  for  0-C0  together 
with  an  indirect  result  from  combination  of  Ar-CO,  Ar-Ar  and  U-Ar 
assuming  that  CO  is  spherically  symmetric. 
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14.  Exptl. 


and  theoretical 


results  for  the  He-Xe,  Xe-He  potential 


Reduced  coordinates  defining  the  interaction  of  an  atom  with  a 
homonuclear  diatomic  molecule  whose  equilibrium  bond  length  is  6/a. 
The  center  of  mass  and  internuclear  separations  are  respectively 
R = p/a  and  r.  = p.  / a. 
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